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ABSTRACT 


The basic governing equations for the study of 
problems of equilibrium and free vibration of shallow, 
isotropic, spherical shells spinning about their polar 
axis are derived from a variational principle. These 
equations are specialized to the case of a spinning flat 
membrane disk and the problem of transverse vibrations 
is solved for the case of a fully clamped hub and the 
case of a loosely clamped hub on an annular disk. The 
method of approach for intermediate hub configurations 
is also discussed. 

The vibration analysis of shallow shells is approached 
and it is noted that the spinning shell equilibrium 
problem must be solved first to provide the necessary 
stress and displacement distributions for the vibration 
equations of motion. The equilibrium problem is treated 
by linear theory for a freely spinning shell and by both 
linear and nonlinear theory for the spinning shell with 
a fully clamped hub. It is found that bending effects 
are important, particularly for shells with a small 
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ratio of rise to thickness, and that a continuous 
transition from the shell solution to the solution for 
a spinning flat disk requires that bending be included. 
The nonlinear effect of finite rotations of shell 
elements is found to be strong, and the effect increases 
as the ratio of shell rise to shell thickness increases. 
It is concluded that a nonlinear theory which includes 
both membrane and bending effects should be employed for 
the general analysis of the equilibrium configuration. 

The equations of motion for small transverse 
vibrations about the equilibrium configuration are 
formulated and reduced to a form suitable for numerical 
solution by an extension of a known method. 
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CHAPTER I 


INTRODUCTION 

The purpose of this investigation is to formulate the 
equations of static and dynamic equilibrium for shallow 
spherical shells which are spinning about their polar axis, 
and to obtain solutions for the static stress distributions 
and free vibration natural frequencies and mode shapes for 
several cases of shell curvature and hub configuration. The 
shells considered in this analysis are assumed to conform 
to the usual restrictions of shallow shell theory and in 
addition are assumed to be isotropic. The static stress 
analysis is formulated allowing for finite rotations, but 
solutions are obtained both with this assumption and with 
the assumption of infinitesimal rotations. The free vibra- 
tion problem is formulated as one of small perturbations 
about the static equilibrium configuration. Various aspects 
of this general problem area have been treated in the 
literature, however the unified formulation of the theory 
and the solutions for the special cases presented here are 
new. 

The necessity for the design of efficient, light- 
weight structures for space applications has led to the 
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consideration of centrifugally stabilized shells. In this 
type of structure advantage is taken of the fact that the 
static stress distribution in a shell spinning about its 
polar axis essentially generates a pre-stressed configuration 
which has a substantial resistence to deformation from the 
equilibrium shape even though the material has practically 
no inherent structural rigidity. In this way a shallow 
shell, such as a radar dish, can be constructed of a flexible, 
light-weight synthetic material and stabilized in shape by 
spinning it about its axis of symmetry. 

Various applications of centrifugally stabilized shell 
structures have been proposed [1,2, 3, 4]. Included among 
these are collectors for solar furnaces, radio and radar 
antennas, solar sails, and drag modulation devices for re- 
entry from space. All these proposed uses have in common 
the necessity to provide a precise shell or disk geometry 
in the presence of known or random disturbances with a 
minimum of structural weight. 

Initial interest in the solution of the vibration problems 

of spinning disks, the limiting case of a spinning shallow 

spherical shell, arose in connection with the design of 

high speed rotating machine elements such as saw blades 

and turbine disks. In general, a complete analysis of these 

problems would require the consideration of the bending 

stiffness of the disk material and the in-plane membrane 

stresses due to centrifugal body forces. Early investigators 
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found that the solution of the general case posed consider- 
able mathematical complexity and elected to study separately 
the cases in which rotational effects were negligible and 
the cases in which bending stiffness effects were over- 
shadowed by rotational effects. Analysis of the situations 
in which bending stiffness is predominant led to a mathemat- 
ical formulation for plate vibrations which was studied 
initially by Kirchoff and which is presently developed to a 
high degree. 

The case in which rotational effects are of primary 
importance led directly to the study of spinning membrane 
disks. Lamb and Southwell [5,6] discuss the problem of 
small transverse vibrations of a spinning circular membrane 
which has no hole at the center and has a free outer edge. 

This problem is shown to reduce to that of obtaining solutions 
to the hypergeometric equation which are finite at two 
singularities. The solutions are in the form of Jacobi 
polynomials. Southwell [6] extends the analysis to the 
study of the transverse vibrations of a spinning membrane 
whose transverse deflections are constrained to be zero in 
an interval r < a and whose outer edge is free. The central 
clamping is such that radial displacements are unrestrained. 

The solution to this problem necessitates finding hypergeometric 
functions which vanish at r = a and are finite at the outer 
edge. Simmonds [7] considers the same problem and points 
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out that Southwell ' s paper contains an error since it over- 
looks the fact that the hypergeometric equation cannot have 
two independent solutions involving only terms containing 
the hypergeometric series 2 F 1 (a, 6, y, x) if the parameter 
Y assumes integral values. Johnson [8] has also published 
solutions to the same problem as part of his investigations 
of the vibrations of shallow elastic membrane shells. His 
results reflect the same error as Southwell's for non- 
axially symmetric vibrations. Lamb and Southwell [5] use 
their membrane analysis results, Kirchoff's non-rotating 
results and Rayleigh's principle to approximate from below 
and above the natural frequencies of free vibrations for a 
spinning disk including both bending stiffness and rotational 
effects. 

The case when the central clamping is such that radial 
displacement is completely or partially restrained is more 
difficult to analyze. The restraint on radial displacement 
gives rise to a more complex in-plane membrane stress 
distribution. The resulting differential equation governing 
the radial dependence of the vibration mode shapes is of 
the Fuchsian type with four regular singular points. In 
principal, this type of differential equation can be reduced 
to a form studied by Heun [9] . However, as discussed in 
the present analysis, it is more convenient not to make the 
reduction in this case. It is found that solutions for the 
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radial mode shapes appear as power series whose coefficients 
depend on a four term recurrence relation. Hence, the 
problem with central clamping is correspondingly more diffi- 
cult in the computational sense than is the problem without 
central clamping. Solutions for the natural frequencies 
of vibration for flat spinning membrane disks with several 
hub configurations constitute a portion of the present 
investigation . 

Bulkeley and Savage [10] have studied the centrally 
clamped membrane for the case of axisymmetric vibrations. 
Their results are for various degrees of partial clamping 
up to and including the fully clamped configuration. 

Simmonds [11] has studied the fully clamped case for axi- 
symmetric vibrations. Both investigations find that while 
the general vibrations are governed by an equation of 
motion which reduces to Heun's equation, the axisymmetric 
case reduces to the hypergeometric differential equation. 

The most difficult boundary value problem for the 
spinning membrane occurs if it is annular and hubless. 

This case was studied by Eversman [12] . It is found in this 
case that reduction of the radial mode shape differential 
equation to Heun's equation is convenient. In this form, 
two of the singularities of Heun's equation correspond to 
the inner and outer edges of the annulus and the boundary 
value problem is that of finding solutions which are finite 
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at two singular points. Eversman also found that the 
axisymmetric vibrations of the annular membrane reduced 
to the hypergeometric equation and that the case of one 
nodal diameter reduced to a form of Heun's equation which 
admits a polynomial solution. This was shown to correspond 
to a rigid body precession mode. 

Interest in static and dynamic problems of spinning 
shells appears to be relatively recent. Explicit attention 
to the equilibrium stress distribution in a shallow spherical 
which is spinning about its polar axis was given minor 
attention by Reissner [13] in his important paper which 
presents solution techniques for the equations of shallow 
spherical shell theory in the case of axisymmetric loading. 

He gives a particular solution of the governing equations in 
the case of the inertia load on a rotating shell. No attempt 
is made to obtain a complete solution for various boundary 
conditions. The shell equations used by Reissner are based 
on the assumption of small displacements and rotations, a 
condition which might not be met if the shell is spinning 
rapidly or the stiffness of the shell material is low. A 
subsequent paper by Reissner [14] based on nonlinear shell 
and membrane theories of Marguerre [15], H. Reissner [16], 
and Bromberg and Stoker [17] allows for the possibility 
of finite rotations. A numerical integration scheme for 
nonlinear two-point boundary value problems for shells of 


- 6 - 



revolution has been presented by Archer [18] . As a limiting 
case of the nonlinear shallow shell theory Simmonds [19] 
obtains solutions for the transverse displacement of a 
normally loaded spinning membrane disk. 

Cohen [20,21] has solved the problem of the static 
equilibrium configuration of a spinning paraboloidal dish. 

His theory is based on a different set of governing linear 
equations than is Reissner's and has the potential short- 
coming that all boundary conditions cannot be satisfied 
simultaneous ly . 

The first investigator who specifically treated the 
statics and dynamics of spinning shallow membrane shells 
appears to have been Johnson [8] . His study is based on 
Reissner's nonlinear shallow shell theory for the membrane 
case. The vibration problem is considered to be one of 
small perturbations about the equilibrium state assumed by 
the spinning membrane shell. He finds that the equilibrium 
state is easily calculated and is one in which there is no 
meridional stress, independent of the curvature of the shell. 
This observation is inconsistent since it is known in the 
limiting case of a shallow shell which is a flat disk that 
the radial stress does not vanish. A consistent theory 
would show a continuous transition from the stress distri- 
bution of a shallow shell to that of a flat disk. This 
continuous transition is accounted for at the expense of 
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a more complicated theory if bending effects are included. 

A second inconsistency in the membrane theory is the 
impossibility of satisfying certain boundary conditions, 
even if the stress distribution is primarily of a membrane 
type. As pointed out by Johnson, this difficulty is probably 
not serious since a boundary layer analysis [14] would show 
that bending effects are important only in the immediate 
vicinity of the edge for shell geometries for which membrane 
theory is nearly correct. 

The general equations for equilibrium of rotationally 
symmetric spinning membrane shells were formulated from a 
variational principal by Simmonds [1]. He obtained the 
equations in a form such that the problem of determining 
the appropriate initial shell shape in the non-spinning state 
in order to obtain a desired shell configuration in the 
spinning state could be solved. As computational examples 
he cites the cases of the flat disk, sphere, and paraboloid 
with conical covering. 

Other than Johnson's [8] investigations of the vibrations 
of shallow membrane shells, there have apparently been no 
published investigations of the free vibration characteristics 
of spinning shallow spherical shells. There have been 
numerous studies dealing with the free vibrations of 
stationary shallow spherical shells. Reissner [22,23] and 
Johnson and Reissner [24] neglect the effect of longitudinal 
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inertia on the predominantly transverse vibrations of 
shallow shells. When this is done, it is found, just as 
in the static case, that the dynamic equations can be 
reduced to two simultaneous equations for the displacement 
and a stress function. Kalnins and Naghdi [25] extend the 
idea to three different degrees of neglect of the effect 
of longitudinal inertia. One of their forms of the equations 
corresponds to those used by Johnson and Reissner. 

The case of vibrations of the spinning shell is 
significantly more complicated than the equivalent problem 
for the non-spinning shell. Of major importance is the 
fact that in the spinning case the equilibrium configuration, 
about which vibrations occur, is not a zero stress state. 

This pre-stressed condition gives rise to membrane restoring 
forces in addition to those associated with the usual shell 
theory. An additional complication arises because of the 
Coriolis acceleration coupling between the components of 
displacement which arises due to spin. Because of these 
additional difficulties most of the solution techniques 
available in the literature cannot be applied to the present 
situation. 

Numerical approaches to the free vibration eigenvalue 
problem for spinning shells are required. Two distinct 
methods appear in the literature. In one method, most 
recently reported by Zarghamee and Robinson [26], values 


for the natural frequency are assumed and the equations of 
motion are integrated by a suitable technique to see if 
the boundary conditions can be satisfied. If they are not 
satisfied, a new value is assumed for the natural frequency 
and the process is repeated. The vanishing of a certain 
determinant which is a continuous function of the assumed 
frequency indicates the satisfaction of the boundary condi- 
tions. An iterative process based on the value of this 
determinant can be used to find the frequency at which it 
vanishes. In the second method, due to Cohen [27], the 
mode shape is iterated on and convergence is obtained when 
the sequence of frequency estimates based on successive 
mode estimates reaches a minimum. The former procedure is 
an extension of Holzer's method and the latter an extension 
of Stodola's method. 

This research program deals with a broad range of 
problems within the general framework of the study of the 
static and dynamic characteristics of shallow spinning 
shells. The general equations of static and dynamic 
equilibrium are formulated from a variational principle. 

The equations of motion for free vibrations are obtained in 
terms of small perturbations about the equilibrium configura- 
tion assumed by the spinning shell . It is shown that the 
linear and nonlinear shallow shell theories of Reissner 
correspond to the present equations in the static case. 
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The spinning disk equations presented by Timoshenko and 
Goodier [28] are the limiting case of these shallow shell 
equations. The spinning membrane shell vibration equations 
of Johnson [8] can be shown to be identical with the 
vibration equations obtained in the present analysis. The 
equations of motion for the transverse vibrations of 
spinning membrane disks studied by Simmonds [7] and 
Eversman [12] and the equations of motion for the in-plane 
vibrations of spinning disks treated by Simmonds [29] are 
special cases of the results obtained in this investigation. 

Several digital computer routines have been written to 
evaluate various aspects of the theory developed herein. 

Three programs have been written to calculate the static 
stress distribution in a spinning shallow shell. One program 
is for the calculation of the stress distribution in a 
freely spinning shell (no hub) using an analytically 
developed solution to Reissner's linear equations. The other 
two programs deal with the case of the shell with a fully 
clamped hub. One utilizes Reissner's linear theory and the 
other utilizes his nonlinear theory. These static stress 
calculations were carried out and are reported here since 
they must be available before a vibration analysis can be 
carried out. Of particular theoretical interest is the 
comparison of the linear and nonlinear theory and the 
shell geometries for which a nonlinear theory must be used 
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and those for which a linear theory is adequate. In addition, 
the shell configurations for which a bending theory is 
required and a membrane theory is inadequate are clearly 
shown . 

The results of Simmonds [7] and Eversman [12] for the 
natural frequencies of free transverse vibrations of spinning 
membrane disks have been extended in this investigation by 
the computation of results for the case of an annular 
membrane with a frictionless hub and for the most physically 
significant case, that of a fully clamped, or "built-in" 
hub configuration. The intermediate case of partial clamping, 
studied by Bulkeley and Savage [10] in the axisymmetric 
case, is also discussed in the general case. 

The equations of motion for the shell vibrations have 
been formulated in a way appropriate for the determination 
of the natural frequencies by the method of Zarghamee and 
Robinson [26]. The procedure is much more complicated than 
their analysis since the equilibrium stress and displacement 
distribution must be computed prior to or during the 
iteration scheme and appears as a variable coefficient in 
the differential equations. No numerical results have been 
obtained due to the limitations of available computational 
equipment . 
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CHAPTER 2 


EQUATIONS OF MOTION 

2.1 Preliminary Considerations 

We will be interested in deriving the equations of 
motion for the free vibrations of a shallow spherical shell 
which is spinning about its axis of symmetry. It will be 
assumed that the shell is isotropic and that the basic 
assumptions inherent in Reissner's linear and nonlinear 
theories [14,30] are valid. The vibrations will be considered 
to be small perturbations about the equilibrium configura- 
tion assumed by the shell when spinning. Since the 
equilibrium configuration must be established before the 
dynamic problem can be solved, the static and dynamic 
problems will be formulated separately. The equilibrium 
problem will be formulated allowing for finite rotations. 

2.2 Axis System and Nomenclature 

The axis system utilized in this analysis is identical 
to the one used by Reissner [30] . As seen in Figure 1, the 
equation describing the shallow spherical shell, measured 
from the base plane, is 
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z = 'H2 z -r x - 


- (.«?- V") 


( 2 . 2 . 1 ) 


and the local slope is 


6 ^. p=r-r _ X 

5c - ^ e 


(2.2.2) 


The approximation in Eq. (2.2.2) is the shallow shell 
approximation based on the assumption that 


R “ 


l 


The basic shell equations will be formulated in terms 
of displacements which are meridional, tangential, and 
normal, as shown in Figure 2. This is the convention used 
by Reissner in his linear theory [30] . In his nonlinear 
theory he uses displacements which are radial, tangential, 
and axial, as shown in Figure 3. The connection between 
the two systems, within the scope of shallow shell theory, 
is 

u> " u> 

H. O. * w> Vfe 

( 2 . 2 . 3 ) 

\J- **. \)* 
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r 


where the barred quantities refer to the system in which 
the deflections are radial, tangential, and axial. 

2.3 Static Equilibrium Equations 

The equilibrium configuration for the spinning shell 
is one of axial symmetry since the shell and load are both 
axially symmetric. The following fundamental relations, 
specialized for the axisymmetric case and consistent with 
the theory of shallow isotropic shells [14,30,42], will be 
used: 

Strain-Displacement 


<L f = 

' f * g * ‘T 

Stress-Strain 
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Curvature Change-Displacement 


V?® = - 7? o3 r (2.3.3) 

The co-ordinate ^ which appears in the strain-displacement 
relations is the distance from the mid-surface of the shell 
to a shell element, measured outward along the local normal. 
The strain energy of the deformed shell is given by 


U (2.3.4) 

If we write the strain-displacement expressions as 

£(■ - e r + 

e e + ^14, (2-3.5) 

where e r and e # are mid-surface strains, we obtain for the 
strain energy 
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tl " Y ' ve © + ^-^e f eel'TcifcieclC 

•v'^vr ^^JjrcJfcieciC 

•+ + ^ (^©* ® o)^r\J Cdrc^Qci 5 ? (2.3.6) 

If we integrate on ^ and note that e r , e^, < r , k , being 
mid-surface strains and curvature changes, are not 
functions of^, we obtain 

U + + •* '?e,e e \ -rdrd© 

+ % I\lt£+v£ + ^>?tf t lQcdl'<Ae (2 - 3 - 7 > 

where 

e.£ 

d -- r?cr3*> 

and t is the shell thickness, assumed constant. 

As shown in Figure 4, the inertia load of an element of 
the shell is horizontal and given by 
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where m is the constant mass density of the shell and w is 
the spin angular velocity. The corresponding potential 
for the element is 

cW-j. - - ^ w\ uO^ r +• U. -V -^_\f df Q 

By combining the strain energy of the shell and the potential 
energy of the load y the total potential energy of the shell 
is given by 

\j ’ ^ l + e 5 £r e ©l 

■v x\y<- T ^ 

- V drde (2 ' 3 ‘ 8) 

By introducing the strain-displacement relations and using 
the principle of minimum potential energy we obtain the 
variational problem 
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& BMW + ^ + k + Wi 2 

+T-A.V* + f + fl] + t^rr 

-V (4: U^A- "Z.\? L ^~^ r 7 - w\ u^(jT*V a+ u^^I^cAq ■= o 

J J (2.3.9) 

The Euler equations for this variational problem are [32] 


;jcC - O (2.3.10) 

^ - T3fU <-V) 

^ (5V v ^©V vv\v>D x (r+a.v^^-^ - o (2.3.ii) 

where Qr r and QT^are the mid-surface stresses. The boundary 
conditions to be satisfied will consist of the hub conditions 
on some hub radius and the natural boundary conditions on 
the free outer edge, which are 
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0> cr + \p ^ - <=> 

-- ° 

S$- •=. o (2.3.12) 

These conditions correspond respectively to the vanishing of 
the bending moment, the Kirchoff condition, and the vanishing 
of the meridional stress at the free edge. 

2.4 Reduction of the General Equations to Special Cases 
A. The Spinning Flat Disk 

In this case Eq. (2.3.10) is the only one applicable 
since when R becomes infinite there is no transverse load. 
Hence as S?-»-«©we obtain 


-v wa (c -v tO) - <=> 


(2.4.1) 


which corresponds to the result obtained by Simmonds [29] . 
In the case when we assume 


^ \ 


we obtain 
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o 


(2.4.2) 


which corresponds to the governing equation derived by 
Timoshenko and Goodier [28]. Simmond's equations for the 
finite deflection of a normally loaded spinning membrane 
can be obtained from Eqs . (2.4.2) and (2.3.10) by allowing 

t^-wO^D-^oand replacing the inertia load in Eq. (2.3.11) 
by a normal pressure -p to obtain 


dr vnrxuj^r^ = O 

i- d / N 

C dr ) ~ i 3 (2.4.3) 

B. The Linear Reissner Equations 

The linear Reissner equations in the axisymmetric 
case can be obtained by neglecting the non-linear term in 
Eq. (2.3.11) which arises because of the inclusion of the 
effects of finite rotations, and by assuming 

7( u. + ^ v 
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With these conditions we obtain 


3c - ST© + r'nuD'V 2 - ' ■=> (2.4.4) 

■q Z V Z C& + "{5 (}SV * Ve) - -O (2.4.5) 

Equation (2.4.5) corresponds to Reissner's Eq. (16) in the 
axisyminetric case [30], with 




yv\ uJ z -r ' Z ~ 


Equation (2.4.4) corresponds to Reissner's Equation (3a) 
in the axisymmetric case with 

~ - >rr\ uO"^ r 


C. Reissner's Nonlinear Equations 

Reissner's shallow shell equations which include 
the effects of finite rotations [14] can be obtained from 
the present theory by introducing a stress function and 
deformation variable, and properly accounting for the 
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transition from displacements in the meridional, tangential, 
and normal directions, to displacements in the radial, 
tangential, and axial directions. 

Prior to obtaining the equations of equilibrium, the 
appropriate strain-displacement relations for the pertinent 
type of displacements can be obtained from Eqs. (2.3.1) 
and the transformations of Eqs. (2.2.3): 

- U r - “iz u$y. -v C l£ r 

^G (2.4.6) 

Here we have not used the barred quantities to distinguish 
the co-ordinate system. 

Equation (2.3.10), with the assumption 


VA, "V \ 


is equivalent to Reissner's Eq. (2.21). In the present 
notation we have 



■V v>r\ u 


o 
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This implies the existence of a stress function such that 


f* - Vsp 

~ -V wuO Z f' z - 

The differential equation which ^satisfies can be 
obtained from the appropriate compatibility equation. It 
is not difficult to verify that we must have for the mid- 
surface strains [14] 

C 4- S’© - ^ - O 

By making use of the stress-strain relations 

e f -- ^ <r©] 

e® -- -k\_^© - >? 

and the definition of the stress function we obtain 

*^\jT HVc * HV “ ° 

If we introduce a displacement variable such that 

d u&- ^ 
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there results 


_L\ in , _ HV“1 x.~ ^ _ (>^ w%uiV 

El^f ^7 - r c - ” - — ■£ 

(2.4.8) 

A second equation is obtained from Eq. (2.3.11) by 
introducing the stress function and displacement variable 
and employing the approximation of Eq. (2.4.7): 

T \ 4vc 4 4>c " ** - ° 

Upon integration we obtain the second governing equation 

*? - » 

(2.4.9) 

where the constant of integration has been shown to be zero 
because of the free edge conditions. 

It is of interest to note that Eqs. (2.4.8) and (2.4.9) 
require only four boundary conditions. The vertical edge 
reaction at the outer edge is automatically satisfied and 
no condition on u^-is required unless the displacement is 
desired from a further integration of<|). 

A second form of Eqs. (2.4.8) and (2.4.9) in terms of 
stress resultants rather than stresses is 
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4d ■*„ » - -t5^r 

- = 

where 

r t <r r - r M rr •* v 

*■ ^©<s> » “f|p -v 

and 

p - p«? v>ki\-t A5?EA = rv-vt 

2.5 Equations of Motion for Small Free Vibrations About 
the Equilibrium Configuration 

The equations of motion for the free vibrations of the 
spinning shallow spherical shell will be obtained by con- 
sidering small perturbations about the equilibrium 
configuration assumed by the spinning shell. The dynamic 
theory will be based on the same assumptions as the static 
theory with the exceptions that the perturbations will 
have infinitesimal rotations and the possibility of asymmetric 
motions will be included. 

The fundamental relations to be used are the following: 
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Strain-Displacement 


£ - u. <■• + ^ ^ ^ ~ + C l^r 

“ $ * 7? * ^ *4( T e T + <’'<'«- 

*ce~ + <* ' $ + ^ e + OC f6 = ^ f6> * <k! re 

( 2 . 5 . 1 ) 

Stress-Strain 

< 3 > -- * ' ?6 ®1 

+ '* >6 <'l 

%© * ^ ^ r ® 


Curvature Change-Displacement 


K,r ~ - uS'v'r 



The strain energy of the deformed shell is 


U * MljW e f + ^e £ © +r retfc®\' rc ^ecJ? 

In terms of strains alone this becomes 

u - k S J[ dr d*)? 

In the static case the inertia load entered the 
formulation as an equivalent applied load and hence as 
a contributor to the total shell potential energy. In the 
dynamic case the inertia load appears naturally in the 
kinetic energy and should not be considered separately. By 
following the development in the static equilibrium case 
we obtain for the total potential energy 

V * \\ C e f T + ej + ^>?e r e e + ^ ^ 

+ + "r^- t ^G©) * T'&r ^T*' 2 - 0 ^©©) 

+ Ct ^ 1 ^drcie 


( 2 . 5 . 4 ) 
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where e . e_, g__are the strains in the shell mid-surface. 

r e r« 

To obtain the kinetic energy we must derive the 
expressions for the absolute velocity components of a shell 
element. In an axis system fixed to the shell and oriented 
with axes in the local meridional, tangential, and normal 
directions the velocity of a shell element relative to the 
shell is 

V* - u e«p + \y e f 

where e^, e^, e r are unit vectors in the meridional, 
tangential, and radial directions, respectively. The 
angular velocity, in the same axis system is 

2o - - 

where the shallow shell approximation has been used in 
resolving the components. The position vector of a dis- 
placed mass element is 

£> " U. + -0- € e +* CA> ef,r* 

and the velocity of the origin of the axis system (coincides 
with the undeformed element) is 

-- r e e 
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The absolute velocity of an element is 


v -- v 0 + 4 cot’s 


which yields 


V " (.& -0O-O-) 4\_o* + (f + a+ OS'—") u>3 e© 

4 ( c> - i>~ oj) e r 

( 2 . 5 . 5 ) 


The kinetic energy of the deforming shell is then written 


T - ” w\^u. - *(C4 0.4 


+ 0 >- " O-^rco'] 


e> 


( 2 . 5 . 6 ) 


The governing equations of motion are obtained from 


Hamilton's principle: 


By replacing e r , e^, and g r ^ in Eq. (2.5.4) by their 
definitions from Eqs. (2.5.1), we obtain for the variational 
problem 








-^*^-¥‘■^• 1 ^- ?{-* 


z 

fr 


+ (io>+ 7-zU^ e ')' Z + -Z'? uv(t + 7* tfi ee) 

+ t>>ejl Z ^3 'f'dfded't m o 


(2.5.7) 
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The shell displacements will be taken in the form 


u 0' , , e i-0 = u e (r) + 

0-d,e,0 - V-*( r,ej) 

- i& 0 (r) + o>%r iQ) i) 

Where cf Jr). Air) are the shell displacements in the 
equilibrium conf iguration of spin about the shell axis and 
u* , v* f o3-* are assumed to be small perturbations about the 
equilibrium position. If these expressions for u, v, and 
U> are inserted into Eq. (2.5.7) and the result is 
rearranged slightly to group terms involving the equilibrium 
displacements, we obtain 

s liliW. + ij '* i 

. tr * * t * m ] 

- iBe.^ 

+-2^^111^^ t3-*-V (_ U*-V oo] J 

- PF-jl u/+ C H ^/] 
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w\> * %Ic * #4 a>; z ] 


~t{ o^tS-^ + 7 uV(r cfy* + 7V.t^) 


+ V 1 oSiy f (7 «^r + 7 l “^ee ) + £ fcWe cH:| 

+ mSSHt^ ■ °-*£*T 4 **-«• «T 4<>M^ «*$-£} 

+ uV^+^T*^* 7 + #4(^)T 

+xv“\_a? + <&>r*Sf +%# 7 ' 7^)4 J 



US £* , 

— - r + 


^V.^)f} - f f «>,** 


* L X r <*? + 1 * &eif *T.J A 4 Cr&r*+ ~k <&J) 

* rjcjejj- o 
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By referring to Eq. (2.3.9) r which defines the 
equilibrium configuration, and by using the fact that the 
processes of variation and integration can be interchanged, 
it is noted that the variation of the first integral in 
Eq. (2.5.8) vanishes. To interpret the second integral 
term in the variation we note that the strain-displacement 
relations of Eqs . (2.5.1) become 




Uo r , . _ . 

<« • t# - v- f- -«•*(+ «v ) 


We note that the strains can be broken down into the sum 
of equilibrium strains and perturbation strains 

+ €* 

~ €-e 0 + 

tir& ~ tire 

where, as previously defined, the equilibrium strains are 

£f 0 " + Z£ U&Of - 7 ^ > L&Off' 

/ J^o I L&b ^ V&o r 

~ r + p? “ ^ “pr 

The perturbation strains, defined as the additional strains 
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due to perturbation of the shell from its equilibrium 
configuration, are 



U. 



(fro f uP* * ^ u^t - <" uP f * 




dl 

r 



--zc 


% U 3 6 B 

hti"*S) 


The second integral term can then be written 

£ JJJ \ ^^00^+ Uo^ u>o£][ <>*+( U*+ 

■*■ ~f + T^^ee +* V 3 r dr ol© d"i 


This is recognized as a variation problem involving the 
equilibrium loads and stresses acting on the perturbation 
strains and displacements. As such, the variation must 
vanish. Hence, the variation problem reduces to 
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The Euler equations for this variational problem are [32] 


0 - <$•© + ~^r§ G - -*z.to i>*- (u*+ ei>*£)co z 1 

(2.5.9) 

+ ^r*S " ri^r — O^u/ 2 "! 

(2.5.10) 
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?vVu/ -*7-3? { r + f <Cf o 

'r& 0 u$g + c&or%e J 

~ - W\^ u>* -ZcO •<>*■£ - ( U*+ UJS^) u/£7 

K ^ J (2.5.11) 

whereQT*, * 3 ^*, 0^* are mid-surface perturbation stresses and 

where, in Eq. (2.5.10), the shallow shell approximation 

r* 

I , was used. The boundary conditions for this set of 
partial differential equations will consist of the hub 
conditions at the hub radius and the free edge conditions 
at the outer edge of the disk. The dynamic boundary condi- 
tions are the same as static boundary conditions and are 
f - 3 - 

u[a) ~o 
O'(a') =o 
u^(a') = o 
u3?(a) =-o 

r - b 

^rr + ^ (~f **■ 7* * ° 

- o 

- o 
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where <**a is the hub radius and r* bis the shell outer radius. 
In terms of the perturbation quantities these become: 

r - a - . 

a* (a} ~ 0*(a) - u9*(a) * o 
uy r *(a) - o 

r - Vd *. 

+ )l r - ^ 

-/z J 1 nV jL /.!• 

v cO* 4- r T>rt ^ "dG'zJlr- b 

<00 - t r *(b) = o 

It should be noted that the displacements are referred 
to the unstressed configuration and hence the perturbation 
quantities are in directions established by the shape of 
the shell before it is spun. Simmonds [1] has circumvented 
this problem to some degree by establishing an intermediate 
reference frame associated with the shell equilibrium 
configuration and perturbations can be referred to this 
system. For vibration problems of shallow shells it is not 
felt that this refinement is justified. 

2.6 Reduction of the Differential Equations for Special 

Cases 

A. The Transverse Vibrations of a Flat Spinning 


- O 


- o 
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Circular Plate 


In this case the radius of curvature is infinite 
and the equilibrium deflection in the normal, or transverse, 
direction vanishes. Equation (2.5.11) decouples from 
Eqs . (2.5.9) and (2.5.10) and becomes 

( 2 . 6 . 1 ) 


The boundary conditions are 

- o 

=o 

u + ^ C f + 7^ ) I ~ ° 

[f? (vM+ ^ ItCt )]| r . o 

B. In-Plane Vibrations of a Flat Spinning Circular Plate 
As noted in (A) , in the case of a flat plate the 
transverse and in-plane vibrations decouple. In this case, 
with the radius of curvature infinite, the in-plane vibrations 
are governed by Eqs. (2.5.9) and (2.5.10) 


- 39 - 



^r(f<5- r *) * vY»r\. u * -•£ olM>* - u*oO*] 


( 2 . 6 . 2 ) 


^O^e) + -|? e + / Ty© - ^<'l^*+'ZcOu' ¥ - <^*<-0**] 


(2.6.3) 


with boundary conditions in the clamped hub case 

M*(s) - o 
<J*(a ) - o 
<3-*(la) -o 
'Tre(b) 

This problem corresponds to the one studied by Simmonds 

[29] . 

C. The Transverse Vibrations of a Flat Spinning 
Membrane Disk 

In the case of a membrane disk the bending stiffness 
vanishes, so that D/t =0. We will thus have 

T^( r< v 0 u>^) 

(2.6.4) 
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which is the governing differential equation for the studies 
of several authors [5,6,7,10,11,12]. 

In this case the boundary conditions are different than 
for the case of the plate because of the vanishing of the 
bending stiffness, and hence the bending moment. At the hub 
radius only the transverse displacement can be set to zero, 
while at the outer edge only the requirement of the finiteness 
of the deflection can be imposed. Hence we will have 

i&Ca) - o 
o>( b) - fi Mite 

D. The Vibrations of a Spinning Shallow Spherical 
Membrane Shell 

In the membrane shell theory, the bending stiffness 
vanishes so that Eqs . (2.5.9), (2.5.10), and (2.5.11) become 

-I ?6V) 

(2.6.5) 

+ ^ £*rt,o(u*+ u>%) -d-V?] 

( 2 . 6 . 6 ) 
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■Tvr { r cr r * + r '3r 0 +T vb (F ' S e« s w ^ + "%r^e J 

+ t ( <5>* 4 ^0 “ cS^- Z aJ> 0 ~*t 9 - (u*+uJ^) CO 7 ^ J 

(2.6.7) 

These equations can be shown to be equivalent to those 
derived by Johnson [8] and solved for some special cases. 

This equivalence can be established by noting that Johnson 
utilizes displacements which are axial, radial, and 
tangential, whereas the present analysis uses normal, 
meridional, and tangential displacements. Within the shallow 
shell approximation the relationship between the two 
displacement systems is 

lA - lA -l- 

0 * - 

uP - uP - uP 

( 2 . 6 . 8 ) 

Furthermore, to the same order of approximation, the stresses 
and are considered to be representative of either the 

mid-surface or "horizontal" stresses. 

It can be seen that Eq. (2.6.6) becomes 
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r >?e ) + ^5® + 'ht = n* - <> V 1 ] 


(2.6.9) 


Equation (2.6.5) can be used to show that 




( 2 . 6 . 10 ) 


By using Eq. (2.6.9), Eq. (2.6.7) can be written 


r*Af (~ ^+^&r)^r f + rSr 0 +7 'aej(~ £ + + r^J^el 


— w^\ 


\ *v* r 7 

o5- - u e ] 


(2.6 . 11) 


The shallow shell relations of Eqs. (2.6.8) can be 
employed to yield 

^ + r<Tr u3- r ] + f? * ^0f) ^T r ^0 o (> % J 


- /vn u&'*~ 


( 2 . 6 . 12 ) 


An equation for dynamic equilibrium in the horizontal 
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direction is obtained by resolving Eqs. (2.6.5) and (2.6.7) 
in the horizontal direction and utilizing the shallow shell 
approximation that < 3 ^ and 'tyg are taken as the horizontal 
stresses. Following this procedure we obtain 


J_ 3L 

r '©r 


) ~f e + r u3- S’ ) 


-“203 t>* - (u* + U&'Q ) 


(2.6.13) 


where it has been assumed that 


r 




Again, employing Eqs. (2.6.8) we obtain 

■=©?( - vys r [ &[ **- "Zco - u v to 2 ] 

(2.6.14) 

In summary, for the case of displacements in the axial, 
radial, and tangential directions, we have, in the membrane 
case 

•ir(i “'Sf) « nnr\a* -"240 - aVj 

(2.6.15) 
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+ + +*3cO a*- t>-*60 z ] 

(2.6.16) 





^T^'^^rW* ^’Te,^] = ws^ 


(2.6.17) 


Equations (2.6.15), (2.6.16), and (2.6.17) are identical 

to Johnson's equations specialized to the case of a spherical 
shell. 

To supplement Eqs. (2.6.15), (2.6.16) and (2.6.17) it 
is necessary to give the mid-surface strain-displacement 
relations for the perturbation quantities. We previously 
have stated these relations for the axis system which is 
meridional, tangential, and normal. They were 


e* - u* + ig* + uhfoS-f <A? Z 


6* - 



tJ* j_/’u5feV 
<? + tl r j 





If we make the transition of co-ordinates given by 
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Eqs. (2.6.3) these become in the radial, tangential, and 
axial directions 



u*+C o-V- 5 ) 4 


X 

z 




In the case that bending is considered we can obtain 
a set of equations of motion for the radial, tangential, 
and axial axis system by including the bending term in 
Eq. (2.6.17) since within the scope of shallow shell theory 
it is unaltered. Hence, in addition to Eqs. (2.6.15) and 
(2.6.16) we will have 


T vV* 6 &*--rlr{ r ("e + ^)'5r* + J 

+ <■&>,•) r re 4 r = - wi 5 * 

(2.6.19) 

The vibration equations in the form of Eqs. (2.6.15), 
(2.6.16), and (2.6.19) have the advantage that the orienta- 
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tion of the axis system is not tied to the initial shell 
shape. 

2.7 Equations for Computation of Bending Moments 

It has already been noted that the vanishing of the 
bending moment arises as one of the natural boundary 

conditions at the free edge of the spinning shell. However, 
it has been felt appropriate to give here expressions for 
the bending and twisting moments in the shell since they 
do not arise naturally in the course of variational 
development. 

By appropriate integrations of the strain displacement 
relations of Eqs . (2.5.1) we obtain the moment resultants 

(moment per unit length) 


tVffY - - of 

7 “V + 'r za %© + 

( i - vO O C ? 

(2.7.1) 

In the case of axisymmetric deformations these reduce to 


- 47 - 



nil I ■■■ III 


M rr - -o^w3* fr + ^ u9fj 

^9© ” " T r + J <A~rJ 

M fe - o 


( 2 . 7 . 2 ) 
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CHAPTER 3 


THE TRANSVERSE VIBRATIONS OF A SPINNING 
MEMBRANE DISK WITH A FULLY CLAMPED HUB 

3 . 1 Introduction 

As pointed out in Chapter 2, the general equations for 
the free vibrations of spinning shallow spherical shells 
contain as a special case, when the shell curvature is 
zero, the equations of motion for the vibrations of a flat 
membrane disk. It is found that the transverse vibrations 
are uncoupled from the in-plane vibrations. The case of in- 
plane vibrations was studied in detail by Simmonds [29] . 

Lamb and Southwell [5], Southwell [6], Simmonds [7, 11], 
Johnson [8], Bulkeley and Savage [10], and Eversman [12], 
have investigated the transverse vibrations of spinning 
membrane disks for several hub configurations and symmetry 
conditions. None of these authors has presented results 
for the case of general vibrations of a disk which has a hub 
configuration of significance for applications, as discussed 
in Chapter 1. It is the purpose of the next two chapters 
to study the transverse vibrations of spinning membrane 
disks which have hubs which provide various degrees of 
central clamping. Two cases are studied in detail. The 
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case to be considered in this chapter corresponds to the 
case of a fully clamped, or built in, hub which constrains 
both vertical and radial displacements of the disk. The 
case to be considered in Chapter 4 is that of a "loosely 
clamped" hub which prevents vertical displacement but 
allows radial displacement. These two hub configurations 
are the extreme cases of physically significant hub 
possibilities. The intermediate hub configurations have 
been studied for axisymmetric vibrations by Bulkeley and 
Savage [10] and their results are used to formulate the 
equations for these hub conditions for the asymmetric case 
at the end of Chapter 4. The hub configurations are shown 
in Figure 5 . 

3.2 The Equilibrium Stress Distribution 

The hub-disk configuration to be analyzed is shown in 
Figure 6. The differential equation which governs the 
equilibrium stress distribution in a spinning disk was 
obtained as a special case of the spinning shallow spherical 
shell. It is given by Eq. (2.4.2) and is repeated here for 
convenience 




- + v'V-i 


~ O 


(3.2.1) 


- 50 - 



In the case of the fully clamped hub the boundary conditions 
to be satisfied are 


u(a) * o 

^(b*) - o 

where r is the hub radius andf*bis the disk radius. 

Solutions for the equilibrium displacement and stress 
distributions are in the form [28] 

U - C, f wvuiV 3 

(3.2.2) 

^ 0 + flc, + 6 -*)% - 

(3.2.3) 

(3.2.4) 

The boundary conditions lead to two simultaneous equations 
for the constants and C 2 
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C,+ 


Qz 

of* 


O-v*) 

«S>E 




(3.2.5) 


b c <3,H 

(3.2.6) 

By solving Eqs. (3.2.5) and (3.2.6) we find for C 1 and C 2 






(3.2.7) 


^ 3 b UiW)^+0-v0a x J 


(3.2.8) 

For the equilibrium stresses we thus have 


(3.2.9) 
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(3.2.10) 



1 /> 4 - 7 

^4 j 


where 


, . (Wj vr, w WL z [ 7 
•""« ' «< LO+vOb^+O 






v\o oJ 


(3.2.12) 




+ ^ V 3 


r^i tO 


(3.2.13) 


3.3 The Governing Differential Equation for Free Vibrations 
The partial differential equation which governs the 
transverse vibrations of a spinning membrane disk is given 
by Eq. (2.6.4) and is repeated here 


r -c>r t ^ V ^ r J + c^&^b ~ ^-£*2- ~ ° 



Where and trr are the equilibrium stress distributions 

given by Eqs. (3.2.9) and (3.2.10). The boundary conditions 
to be satisfied for the fully clamped hub are 

c a-£a) - o 

ca-Lk,) - f.n i-re - 

Hence we seek solutions of Eq. (3.3.1) which vanish at the 
hub and are finite at the outer edge. 

We assume a solution of the form 

Qji) - vV(Y) -sin 

which, when substituted in Eq. (3.3.1), yields the following 
ordinary differential equation for W(r) 


X <£. 
Coir 




- CD 


(3.3.2) 


where 



xi 


It is convenient to partially non-dimensionalize Eq. (3.3.2) 
by introducing the non-dimensional radius P defined by 
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r ^ pb 


If we introduce this into Eqs. (3.2.9) and (3.2.10) for 
the stresses we obtain 


(3.3.3) 



(3.3.4) 


where 


•^-VV J L b/ L^I 4- V^)lcT^ v(.| 

(3.3.5) 


With the introduction of the non-dimensional radius and 
Eqs. (3.3.3) and (3.3.4) into Eq. (3.3.2) there results 

j$\j 0 f'*r'^r$h 

(3.3.6) 
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Further simplification can be obtained by introducing the 
change of variable 



and the definition 




1+lylz , 3z. \l^± Ojjvj 

■ 54 V ^ •%Vv'U.“V S 



to obtain 


KOt-iXx+S*) ^x. +b-X t '-0 -% z )*] "dx + ju/- o 

(3.3.7) 


In this form the boundary conditions to be satisfied are 

wUftfl -o 


W(|) - ITB 


(3.3.8) 


3.4 Solutions of the Differential Equation 

The basic differential equation for the radial dependence 
of the vibration mode shapes is given by Eq. (3.3.7) and 
supplemented by Eqs. (3.3.8) for the boundary conditions. 
Equation (3.3.7) is of the Fuchsian type with four regular 

o 

singular points at x=0, x=l, x=- v and x= oo. It can be 
reduced by an appropriate transformation of the dependent 
variable to Heun's equation [32], In the particular case 
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at hand, direct reduction to Heun's equation, at least 
by conventional means, does not appear practical since it 
involves a transformation of the form 

WOO 

in which k is complex for vibration modes with more than 
one nodal diameter. Another transformation which will 
lead to Heun's equation has been used by Eversman [1.2] in 
the case of the annular membrane with free edges. In this 
approach the transformation 

\ 

** * X 

is made and then a transformation of the form 

w«) * 

This reduction is successful in that Heun's equation can be 
obtained with real values of k. However, the transformed 
geometry of the singular points must be considered. After 
the change of independent variables there will be singular 
points at 1 , , and . The region 

of physical interest lies between * l and < . 

Since the distance between ^ * (V a f and \ is less than 
1 only if Otff < 2 , and since series solutions about ^ * V 

can be guaranteed to converge only for , only 

certain disk geometries can be conveniently treated. 
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Because of these observations and restrictions, no further 
reduction of the basic differential equation is made. 

Solutions to Eq. (3.3.8) can be obtained in the form 
of a power series expansion about one of the singular 
points [33] . The singular point at the free edge is the 
only one about which an expansion can be made which can be 
proved to converge for the entire physical region for any 
disk geometry. Hence it is appropriate to shift this 
point to the origin with the transformation 


(3.4.1) 

where 

\ 

*** * <-(%? 

(3.4.2) 

With this transformation the singular point geometry changes 
as follows : 

V - I »•><-. o 

y o i-' 
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The hub radius, originally at v * t*y , transforms to \ . 
The differential equation transforms to 

-^h 7 - -yU^-dSf] V/ = © 

(3.4.3) 


and the boundary conditions become 


W(o) - ^\i4\-re 

WCO * 0 (3.4.4) 


Solutions to Eq. (3.4.3) are sought in the form 


^CO ' § C„ < 

r\*o 



(3.4.5) 


If we introduce the definitions 
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£ 

+^**0 + *1 + p ) — 

(3„ - - «* Jo+ S, 7 ^ fp^+p- 0 + ( ^ +^X ^ + /°) - i il^*- $*(>-&)] / 

* « 4 \^S t X«+( 9 X* , Y’' l ^^ rt ^X rt+ rt + ?[('- 7 ^)' sr >*' T 7f 
^ «-^(» + S 7 X^^ 

it is found that the unknown coefficients, C n , are defined 
by the recurrence relations 

K c o * ® 

*o C.4- S,c, - o 

+ X.C, + S z C t •* O 

*** ^ ^ yi+l ^ ^ ^ 

The first of these relations is the indicial equation. For 
arbitrary C Q it follows that 


or 


S 


© 


» O 


- qr < ’(i+s t )^ z - o 
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This yields the two indicial exponents of the solution 


P, ' ° 

P*'° 

With |0 s. o there exists a solution 

w-3c. 

y\“o 


( 3 . 4 . 6 ) 


with the C n defined by 


VoC. + s ,c, - 

fif. * < C, ♦ S,C, 


” o 


+ ^ +l C„ 4t - o ; 


( 3 . 4 . 7 ) 


where 



(5^ -- -Of^+S 1 ) "(V1-.) 4(*+&) „ - &*)] f 

- a + \(.W) «(«-.') H4+ZE?) *+i[(i--z?H 7 :x ^l[ 

&,_ - - a fe, Ci+^ t ) vi 7 - 

(3.4,8 

Since the characteristic exponents of the solution are 
repeated, only one solution of the above form exists. The 
second independent solution is known to have a logarithmic 
singularity at ^■*0 [33] and is discarded to satisfy the 
boundary condition at that point. 

3.5 The Eigenvalue Problem 

The boundary condition of finiteness at*£«0 is 
satisfied by the solution given by Egs. (3.4.6), (3.4.7) 

and (3.4.8). The eigenvalue problem is that of finding 
values of the parameter , and hence the natural frequency 
-fO , for which a solution of this form will satisfy the 
boundary condition WO)- o with C Q not zero. By noting 
that C)^ and are functions of we conclude that 
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Hence the eigenvalue equation is 


W(l) ■ X CJyU*) ' O 

v\*o 

( 3 . 5 . 1 ) 


and we seek to find the zeroes of the transcendental function 


H M 1 ) 


2 C^CycA) 


-o 


( 3 . 5 . 2 ) 


3 . 6 Numerical Evaluation of the Eigenvalues 

A digital computer program has been written to find 
the roots of the eigenvalue equation 


FiyU*) » O 

( 3 . 6 . 1 ) 
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The method employed consists of evaluating for a 

sequence of values of until a change in sign is 

noted. When the change of sign occurs, and the presence of 

a root is thus indicated, an iteration procedure based on 

the method of secants [34] is employed to obtain the 

eigenvalue. This procedure is carried out for a given 

number of nodal diameters s and for as many values of S in 

I as specified. For specific values of s and 

(^) any number of roots can be calculated, beginning with 
b 

the smallest and proceeding in ascending order. 

After a particular eigenvalue is obtained, the eigen- 
function can be evaluated, by program option, by evaluating 
the solution for the known eigenvalue , at as many points 
in the interval ©*£”4 I , or lo ,as necessary for 

definition. 

The above numerical technique was programmed for the 
IBM 1620 Digital Computer at Wichita State University. It 
was found that accuracy and convergence of the iteration 
scheme could not be maintained without the use of double 
precision arithmetic. This, coupled with the inherent 
complexity of the evaluation and iteration calculations, led 
to computational times per root which were impractical for 
production computation on the 1620. Through the kind 
cooperation of the IBM representative access was made 
available to an IBM 360, Model 65, as well as a Model 75, 
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for the purpose of completing the numerical work. Computa- 
tion times, on these machines were only a few seconds per 
root . 


3.7 The Special Case of Symmetric Vibrations 

In the special case of symmetric vibrations the 
parameter S-O and Eq. (3.3.7) reduces to 

(jr-.Xjf*?).® - o 

(3.7.1) 

where 


s* 


Z 



(3.7.2) 


If the transformation 

\ - y 

< ’ A + 

is introduced into Eq. (3.7.1) there results the differential 
equation 
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<«™ % -hc-1 0 -&c£fw - o 


(3.7.3) 


with the boundary conditions 


WCcO - 




l 


X - 

\ + V- 1 


O 


(3.7.4) 


The mathematical problem specified by Eqs. (3.7.3) and 
(3.7.4) was studied by Simmonds [11] in his investigations 
of this case. 

It is found in this case that there exists one solution 
of the form 


w(c) 


oo 


v\ ~ o 


(3.7.5) 


where the Q^are defined by 
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r 


’ rx-H 


_ 'T'C V“\ -V 


(3.7.6) 


The second solution is logarithmically singular at ^ •% o 
and is excluded to satisfy the boundary condition at that 
point. 

7 

The eigenvalue problem is that of determining yU such 

that 


. ?cx4^T - « 

V'l'a. o 

(3.7.7) 


so that the second boundary condition is satisfied. As a 
means of verifying the results of the general eigenvalue 
problem specified by Eq. (3.6.1), a digital computer program 
for the special case of symmetric vibrations was developed. 
This program utilized an iteration technique similar to the 
one discussed in Section (3.6) and was based on the eigen- 
value problem specified by Eq. (3.7.7). The results of 
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these calculations showed good correlation with Simmonds 1 
results and agreed perfectly with the results for the case 
^ •» © from the general problem. 


3.8 The Special Case With One Nodal Diameter 

In order to further verify the analysis it was decided 
to also treat independently the case of vibrations with one 
nodal diameter. In this case we have s * \ and the governing 
differential equation for the radial mode shapes is 

xCx-OCx+S*) 5$- +Ltx t -0-?)*] 

-%I fw' =o 

(3.8.1) 

where 

■z j*. \ + '±iy i 


(3.8.2) 


In this case it is convenient to simplify the differential 
equation by introduction of the change of variable 
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(3.8.3) 


W(x) - yc*) 


With this substitution Eq. (3.8.1) becomes 

*(x-»Xx+ ft + {(x-iXx+tf) + x(x-i) + x(x+& z )|^ 

+ ^(^> v^)XY - o 

(3.8.4) 

Equation (3.8.4) is a form of Heun's Differential Equation 
[35] . 

As in previous cases, it is convenient to expand the 
solution about the point X * l . This point is shifted to 
the origin by the transformation 


K 


x- » 

CW - * 


“ 0( z (x- 0 


(3.8.5) 
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The differential equation becomes 


+ ^C\<-O( 7 0+S t )] + <(<-<**) + (<-<**)[< -c^(i+£ T )]^ 

+ - o 

(3.8.6) 


and the corresponding boundary conditions, deduced from 
Eqs . (3.4.4), are 


Y (.o) - m v-rc. 

Yco •- O 


(3.8.7) 


One solution to Eq. (3.8.6) exists in the form 



(3.8.8) 
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where the recurrence relation for the C is 

n 


+ (3 0 C 0 •=> o 


(3.8.9) 


and 


fi ^ - - c* z lU +?) •v'C.w-O} 4 (A+'fc? ) y% 4 


If, - 0( 4 (l 4 ^) ^ 


(3.8.10) 


A second solution exists which is logarithmically singular 
at if » o . The boundary condition of finiteness at the 
origin requires that this solution be discarded. 

At this point it is interesting to note the simplifi- 
cation introduced by the transformation to Heun's Equation. 
The recurrence relation given by Eqs. (3.8.9) contains only 
three terms, but it would have contained four terms if the 
transformation had not been accomplished. A similar result 
would have been observed in the general case, but as 
previously noted, other difficulties precluded the transfor 
mation . 
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In this case the eigenvalue problem is one of 
determining jL4 such that 


U Cy*( Z ) * ^ C^) » O 

Y\-<P 


(3.8.11) 


The roots of this transcendental equation were found by a 
method similar to the ones employed for the symmetric and 
general cases. The results were compared with those 
obtained from the general program for the case Sr ) and 
were found to be in complete agreement. 

3.9 Results 

The numerical procedures discussed in the previous 
section were used to calculate the first four eigenvalues, 
v*\®© to vx®^, for cases of from zero to three nodal diameters, 
S-O to . The results of these calculations are plotted 

in Fig. 7 as a function of the ratio of the hub radius to 
the disk radius. 

The general trend is for the increase of jU , and 
hence the vibration frequency, with an increase in the hub 
size. This is reasonable from a physical viewpoint since 
increasing the hub radius would appear to stiffen the 
system. 

An exception to this observation occurs for the higher 
modes corresponding to a given number of nodal circles n. 
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It is observed, for example, that for v\«o and and 

that there is an initial reduction of yli* with increasing 
S . It is anticipated that this trend will appear 
whenever the number of nodal diameters is substantially 
larger than the number of nodal circles. This trend can be 
explained in terms of the dependence of the local disk 
stiffness on the nodal geometry. The equilibrium stresses 
in the disk are given by 


«> - '&o-weo 

r (3.9.D 


(3.9.2) 


where 


c? -ft (v«’) l » 1 -(Hv > )a ' 1 l 

(3.9.3) 


We note that 
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( 3 . 9 . 4 ) 


- v '-f ’ ° 

( 3 . 9 . 5 ) 

so that for the local radial stress increases with 

increasing S while the tangential stress decreases with 
increasing £ . Since & depends directly on (g-) ^ the same 
can be said about increasing or decreasing the ratio of the 
hub radius to the disk radius . 

If there are substantially more nodal diameters than 
nodal circles the local membrane stiffness will depend 
predominantly more on the tangential stress than on the 

radial stress for small (g) and will tend to decrease with 

sl a. 

increasing (g-) up to a limiting value of (g-) . It would be 

expected that the more the number of nodal diameters exceeds 

the number of nodal circles, the more pronounced this 

effect would be and the larger would be the range of (g) 

over which it is observed. 

If, on the other hand, the nodal geometry tends to 

make predominant in governing the local stiffness, then 
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the trend of increasing yrf*’ with increasing (g-) is expected. 

A second interesting feature occurs for (g-) near 
unity where for the the higher values of n there is little 
variation of yt/ with s. This can be explained by physically 
reasoning that for large n and large (g-) the stiffness is 
governed almost completely by the nodal circle geometry since 
they are very closely spaced while the nodal diameters are 
widely spaced in this region. For large s it would be 
expected that more variation with s would be seen. 
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CHAPTER 4 


THE TRANSVERSE VIBRATIONS OF A SPINNING 

ANNULAR MEMBRANE DISK WITH A FRICTIONLESS HUB 

4.1 Introduction 

Treated in this chapter is the problem of the transverse 
vibrations of a spinning annular membrane disk with a 
frictionless or loosely clamped hub. By this it is meant 
that the hub does not restrain in-plane motions of the 
disk but prevents transverse displacements. From a physical 
standpoint this case is similar to the case discussed and 
solved by Simmonds [36], except here the disk is annular 
rather than solid, so that the presence of a hub is on a 
more physically realistic basis. The author was assisted 
by Mr. Thomas Gilley [37] in the following analysis. 

4.2 The Equilibrium Stress Distribution 

The differential equation and boundary conditions 
governing the equilibrium stress distribution in the case of 
a loosely clamped hub are 

Jr 

cry (.s') - o 
C'O - ® 
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(4.2.1) 



where a and b are the inner and outer annulus radii. 

Solutions to the differential equation are given by Eqs. 
(3.2.2), (3.2.3), and (3.2.4). By employing the two boundary 
conditions the two constants, and C 2 , can be obtained 
and the solutions for the stresses are found to be [28] 


<3V - 




(4.2.2) 


* 


WVOO 




(4.2.3) 


In a partly non-dimensional form, with 


r - pb 

Eqs. (4.2.2) and (4.2.3) can be written 


(4.2.4) 


(4.2.5) 
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1 (4.2.6) 

where 

c . d. 

b - b 

(4.2.7) 

is the annulus radius ratio. It is interesting to note 
the similarity between the functional form of the above 
stresses with those given by Eqs. (3.3.3) and (3.3.4), 
keeping in mind, however, the difference in the definition 

of £ . 

4.3 The Governing Differential Equation for Free Vibrations 
By substituting the stress distribution given by 
Eqs. (4.2.6) and (4.2.7) into the governing differential 
equation for transverse vibrations, given by Eq. (2.6.4), 
separating variables, and using the non-dimensional radius 
defined by Eq. (4.2.4), there is obtained for the radial 
mode shapes in the present case 
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+{(^ + ($f+ x )-*'lo+e>p © 

(4.3.1) 

With the change of variable 

y-^ 

(4.3.2) 


Eq. (4.3.1) becomes 

Kfr-iXx-i 1 ) + [ix z -(i+i z )x] 

+i|L(' +P) + £ I s* -^ z x]jw - o 


where 




-34 






W 


(4.3.4) 


The boundary conditions imposed by the hub on the 
transverse vibrations of the disk require 
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W(c) - o 

W ( k ) - P'imite* 


where r*C is the radius of the hub. In terms of the 
variable x these conditions become 


W(-§) = w'Ct*) » o 

(4.3.5) 


WCi) - F'imitE 


(4.3.6) 

where is the square of the ratio of the hub 

radius to the disk radius. 


4.4 Solutions of the Differential Equation 

The observations of Section (3.4) concerning the 
difficulties of reducing the governing differential equation 
to Heun's equation are also pertinent for the present case. 
For this reason, Eq. (4.3.3) will be treated without further 
reduction. 

An expansion in a power series about the singular point 
at X-l is facilitated if we shift this point to the origin 
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with the transformation 


zf - -«*(*- 0 


where 


(4.4.1) 


, _j i 

01 ' i-C*Q z = rj* 

(4.4.2) 

Accordingly, the singular point geometry shifts as follows: 

x - l =*> < * o 

K ’ ^ -=-■» < * Of^O 

X ■*- O '=•'> < = C( z 

£ 

The hub radius, originally at X * £ , shifts to z£« / . By 
introducing this transformation into the differential 
equation we obtain 


+ i [[- a*6 + & T X<-a ? y * ot*l z ]& o 

(4.4.3) 
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r 


The boundary conditions become 

W to ) - Finite- 

WO) - o 

(4.4.4) 

It is found that there exists one solution of the form 

w - t C,*" 

(4.4.5) 

where the C are defined by the recurrence relation 
n 

S 0 C 0 4 E,C, - o 

fr,C«> + y* , c, + S-jC-t = O 

(4.4.6) 

The coefficients fir\ ' ^ are given by 

Of *- ^ ~^u z /a 

- “ W- 
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r 


A second independent solution of the form of Eq. 

(4.4.5) does not exist because the characteristic exponents 
of the solution are repeated. The second independent 
solution is logarithmically singular at j^*oand hence is 
excluded to satisfy the boundary condition of finiteness 
at that point. 

4.5 The Eigenvalue Problem 

The power series solution given by Eq. (4.4.5) and 
supplemented by Eqs. (4.4.6) and (4.4.7) satisfies the 
condition of finiteness at o . The boundary condition 
of no displacement at the hub, | , will be satisfied 

if 

Wo) 3 % C^CyOf 1 ) - O 

r\~o 

(4.5.1) 

Equation (4.5.1), with the defined by the recurrence 

relation in Eqs. (4.4.6) and (4.4.7), constitutes a 

7 

transcendental equation defining the values of for 

which WcO vanishes with C Q not zero. Hence the eigenvalue 
problem is that of finding the zeros of the function 

oo 

vx= ° (4.5.2) 
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A numerical technique similar to the one described in 
Section (3.6) was used to numerically evaluate the roots 
of Eq. (4.5.1). The present case has the added complication 
that in addition to nodal geometry, defined by integer 
values of s, we must consider variations in both the ratio 
of the hub radius to the disk radius , € , and the annulus 
radius ratio, S . 

4.6 The Special Cases of Symmetric Vibrations and Vibrations 

With One Nodal Diameter 

As in the case of the fully clamped hub, there are two 
special cases which reduce the complexity of the governing 
differential equation and can be used independently of the 
general development for computational purposes or for 
verification of results computed from the general case. 

The case of symmetric vibrations, with leads 

to the differential equation 

xfr-iXx-S*) jjji ' ° 

(4.6.1) 

where 



By introducing the change of variable 
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V — * 


< * 


we obtain for the differential equation and boundary condition 





- O 


( 4 . 6 . 2 ) 


V\/C o') - rvMi re 



( 4 . 6 . 3 ) 


It is found that the solution satisfying the condition of 
finiteness at the origin is of the form 


*«) ' § 

■A * © 

( 4 . 6 . 4 ) 


with 


c 

v - y-*4-| 


( yg-v Q - p 

Cvg+/^ 


( 4 . 6 . 5 ) 


The eigenvalue problem is that of determining such that 
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o 


c,(^) - fj cjy)[ jiftJ ' 

( 4 . 6 . 6 ) 


In the case of one nodal diameter, the differential 
equation is 


With the change of dependent variable 

*/ - x^yoo 

and the subsequent change of independent variable 


< » ffji - <x z (i-x) 

the differential equation becomes 


( 4 . 6 . 8 ) 
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<(<- «'•')[ <-d t Ci- S')] 3 f^!* 

+ {<U-o( t o-s t )] + <(<-0^) +(<-a')[c- 

+ i('a>>>'y<-W Z )V “ O (4.6.9) 

and the boundary conditions become 

V(©} - riM»TE 

YCO - o 

(4.6.10) 

The solution of Eq. (4.6.9), a form of Heun's equation, which 
is finite at the origin is 

y<o -3 c R c" 

v\-o (4.6.11) 

with the C n defined by 

AA + <c. = o 

«4-,C«M + AA> + K'y 1+ ,C„ +1 *o; y> V I 

(4.6.12) 

where 
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- >n(v^+Z') + A ( 


Prx ' “ fc-Sft Y\^^\) +(4 -Ztf') ^ -V 

(4.6.13) 

In this case the eigenvalue problem is that of determining 
yiS*- such that 


oo 

UiM*) - "t, C„(ytf z ) - o 

v\~o 

(4.6.14) 

Numerical results for the eigenvalues for these special 
cases were obtained using methods similar to those described 
in Section (3.6). These results were compared with results 
obtained from the general program in the special cases ^■so 
and *3*j and complete agreement was found. 

4.7 Results 

The three Fortran computer programs referred to in this 
chapter were used to compute the first four eigenvalues,^* 0 
to , of the frequency parameter, ytf z , for values of the 
annulus radius ratio, S , from 0.1 to 0.8, for values of 
the hub to disk radius ratio, 6 , from 0.2 to 0.9, and for 

values of s from 0, to 3. The resulting data from these 
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calculations were plotted and are presented in Figures 8 
through 19 . It should be noted that no calculations were 
carried out for cases in which 6 was smaller than or equal 
to & . This was considered a physically unrealizable 

case since it means that the center hole in the membrane 
is larger than or equal to the hub radius. For this reason 
€ was arbitrarily picked to be at least 10 percent 
greater than £> . The limits on £ and S were governed 

somewhat by considerations of practicality in the computer 
programs. For values of S less than 0.3 the convergence 
of the power series was extremely slow and required very 
long computer runs. For values of 6 greater than 0.9, 
the fraction becomes large and leads to possible 

computational difficulties. 

The plotted data. Figures 8 through 19, show that if 
the parameters & , £ , and s are held constant the natural 

frequency of vibration increases as n increases. This 
behavior is predictable from the general theory of Sturm- 
Liouville eigenvalue problems. 

The plots also show that the frequency of vibration is 
increased as the ratio of hub radius to disk radius, £ ; is 
increased and S , s, and n are held constant. The increase 
in the size of the hub with constant nodal geometry and disk 
size causes an overall stiffening of the disk, raising the 
natural frequencies of vibration. 
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In general, the frequency decreases as the annulus 
radius ratio, £ , is increased and s, n, and € are constant. 
The exception for this rule is found in the case where the 
number of nodal diameters is much greater than the number of 
nodal circles. For this case, the frequency increases 
when is increased. This observation was also made and 
explained by Eversman [12] for the case of the annular 
elastic membrane with free edges . 

The phenomenon is shown to be due to the change in the 
local stiffness in the membrane with changes in the nodal 
geometry. The stresses were known to be due only to the 
centripetal accelerations induced by the rotation. The 
stresses are given by Eqs. (4.2.5) and (4.2.6). These 
equations are repeated here for convenience: 



We note that 
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dCS-Q 


') < ° 

w,o=?- <2ge> ^Go*+0 >° 

^ p 


so that for the radial stress decreases with 

increasing S while the tangential stress increases with 
increasing & . 

For the case when the number of nodal diameters is 
substantially greater than the number of nodal circles the 
tangential stresses are predominantly more important in 
providing stiffness and the local stiffness in the membrane 
will increase with increasing 6 . The increased stiffness 

causes increased frequency of vibration. If the number of 
nodal circles is substantially greater than the number of 
nodal diameters the radial stresses are dominant and the 
natural frequency tends to decrease with increasing S [12]. 
If Figures 11, 14, and 17 are studied briefly it will be 
noted that the frequency of vibration is, in fact, influenced 
less by increasing or decreasing the hub radius in the 
cases where the number of nodal diameters is greater than 
the number of nodal circles. This observation can be 
explained. The larger number of nodal diameters causes 
sufficient restraint on the deflection of the membrane in 
thfe proximity of the inner edge that the increase or decrease 
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in the hub size has little effect on the stiffness of the 
membrane or on the frequency of vibration. 

4.8 Intermediate Cases of Central Clamping 

The cases of full clamping and frictionless clamping 
represent the extremes of physically significant hub 
configurations. The fully clamped case corresponds to a 
disk built into the hub or shaft. It can be viewed as 
being equivalent to a disk clamped between collars with 
sufficient clamping pressure to prevent any radial 
deformation of the disk in the region of the hub. On the 
other extreme is the case of frictionless clamping in which 
there is no clamping pressure (alternatively, no hub-disk 
interface friction) and hence no constraint on the radial 
disk deformation in the hub region. Between these two 
configurations are the cases of intermediate clamping in 
which there is sufficient clamping pressure and hub-disk 
interface friction to prevent radial displacements of disk 
elements over a portion of the hub region, but not over the 
entire region. These cases have been studied in detail 
in the axisymmetric case by Bulkeley and Savage [10] and 
the extension to asymmetric vibrations constitutes only 
minor modifications of the procedures previously established 
in this report. 

Bulkeley and Savage show that the form of the stress 
distribution in the spinning disk depends explicitly on the 
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clamping conditions. In particular, for a disk of thickness 
h clamped by a hub of radius c with clamping pressure P 
and Coulomb friction coefficient yl4 r the radius r=g which 
limits the region of zero radial displacement is given by 





with 




For values of there is no radial displacement while for 

<^£1*60 the clamping pressure is insufficient to prevent 
slippage and radial displacement. 

If the radius of effective clamping is between the 
inner radius of the disk and the radius of the hub, 
cJ £ ^ C. , the stress distribution is given by 


r'Xr 1 * c,c*) 


( 4 . 8 . 1 ) 




^ :^ Z ^(V£-€,C r )^-€,cV- 


S + V 3 r J 


( 4 . 8 . 2 ) 


where 
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t *»if 

< * Wi ff 

( 4 . 8 . 

and 

a = inner annulus radius 
b = outer disk radius 
c = hub radius 

g = radius at which clamping pressure and disk-hub 

interface friction become insufficient to prevent 
radial displacement of disk elements 

S'- Coulomb friction at disk-hub interface 

P = clamping pressure 

h = disk thickness 

It is also possible that the effective clamping radius 
will be less than the inner disk radius. If this occurs 
there will be slippage over the entire hub-disk interface. 
This will occur if 

r * 

and the stress distribution will be 


( 4 . 8 . 4 ) 


- 94 - 



<s- e « ^ ^ t*j 

( 4 . 8 . 5 ) 

where 


rtfflrtM 

+4S&4 






( 4 . 8 . 6 ) 


It is recalled that in the case of frictionless 
clamping the stress distribution is given by 


-3-c - 'Hf 

( 4 . 8 . 7 ) 

^ \c ^ ^ - K 4$? 

( 4 . 8 . 8 ) 

where 


& z - (%? 
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Hence/ if the computer program for the loosely clamped case 
is modified to replace 

(4.8.9) 

and the input to compute € is provided, the vibration 
problem for the intermediate case with the effective clamping 
radius less than the inner annulus radius can be analyzed 
by using in Eq. (4.8.9) and the case when the clamping 

radius is between the inner annulus radius and the hub 
radius can be analyzed by using , in Eq. (4.8.9). 

Note that in the case P=0 we find 



Since we replace 

- -(tK - df 

we find that ^and ^^are properly defined for the loosely 
clamped case. 

In the case g=c we have 

c \ ^ - t 1 ±f) C 2 1 

4- O+vO fc? 1 
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and we replace 



_ 1-^/CVtr 1 

’sWK'o) lo-^C^+C^+Ob 7 -] 

By referring to Eqs. (3.2.9) and (3.2.10) and accounting 
for the pertinent definition of for this case given by 
Eq. (3.3.5) [note that the hub radius in this equation 
is a instead of c] we see that and will be properly 
defined for the completely clamped case. 

Hence it is concluded that the complete range of hub 
conditions between, and including, the extreme cases can 
be analyzed by a simple modification of the basic analysis 
for the loosely clamped case. 
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CHAPTER 5 


THE EQUILIBRIUM STRESS AND DISPLACEMENT 
DISTRIBUTION IN A SPINNING SHALLOW 
SPHERICAL SHELL 


5 . 1 Introduction 

In Chapter 2 it was shown that the equations of motion 
for the small transverse vibrations of a spinning shallow 
spherical shell about its equilibrium configuration require 
the knowledge of the equilibrium stresses and displacements. 
It is the purpose of this chapter to present methods for 
the calculation of the equilibrium conditions for the cases 
of a freely spinning shell and a shell with central clamping. 
The linear and nonlinear theories of Reissner [30, 13, 14] 
are employed for this purpose. 

In addition to providing calculation schemes , an 
evaluation is made regarding the adequacy of the use of 
membrane shell theory for the determination of the stresses 
and displacements. This was prompted by the paper of 
Johnson [8] in which was treated the problem of transverse 
vibrations of spinning shallow spherical membrane shells. 

In his analysis it is found that as the curvature of the 
undeformed shell is reduced, that is as the shell approaches 
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a flat plate, the stress and displacement distribution, and 
hence the vibration natural frequencies , do not approach 
those of the flat plate. It is shown in the present analysis 
that this anomaly is due to the membrane assumption and that 
an extremely shallow spinning shell theory (approaching a 
flat plate) must include bending effects if a continuous 
transition from shell results to plate or flat membrane 
results is to be obtained. Furthermore, the inclusion of 
bending effects permits the specification of physically 
significant boundary conditions at the hub and outer edge. 

5.2 The Equilibrium Equations for the Linear Theory 

Equations (2.4.4) and (2.4.5), derived from a variational 
principle, correspond to Reissner's linear results [30] in 
the axisymmetric case with surface loading due to the 
centrifugal force effects. By introducing a stress function 
such that 

^rr *• ~ r 

+ & 

where 

-ft ~ - 4 p<oV* 

and by utilizing the compatibility equation of Eq. (2.4.7) 
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we can obtain two simultaneous ordinary differential 
equations in the deflection and stress function 


2 "Z (\ - V 5 ) 

(5.2.1) 


o ^ 


where 


V 


z 


■5 


£ j. dL 

+ r jr 


and 


(5.2.2) 


(° 


mt = mass per unit surface area 


The boundary conditions for the 
spinning shell are that the stresses 
finite at the point r=0 and that the 
free edge, r=b, are 


case of the freely 
and displacements be 
conditions at the 
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bJ rr - O 

3r (v* "*■ <=> 

d Z oSr , Vi dl& 

37i r dr " 0 

(5.2.3) 

there is a fully clamped central hub of radius a 
add to the free edge boundary conditions of Eqs. 
the hub conditions 

u>( a*) » o 
w ( a) - e> 
u%.(a) « o 

(5.2.4) 

5.3 The Equilibrium Equations for the Nonlinear Theory 
The nonlinear theory of Reissner allows for the 
possibility of finite rotations of shell elements. As 
noted in Eqs. (2.4.6) and (2.4.9), the equilibrium equations 
are 


If 

we must 
(5.2.3) 
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a 


(5.3.1) 


Ji i? 4>V> . tp 

D l«Vr + T~ T + T ’ 

TfJ Wfc. f £•- £f- r+kf'- ^l£^£ 

' 4 6 (5.3.2) 

where 

**rr * ,i ' 5 r * 'r 

* V»* +(><Jf x ' 

?$' + 

The same boundary conditions apply in this case as do 
in the linear case, except that the vertical edge reaction 
(Kirchoff condition) is automatically satisfied and specifi- 
cation of u> is not required unless it is desired to 
obtain uy by integration. 

5.4 Stresses and Displacements in the Freely Spinning Shell- 
Linear Theory 

Because of the lack of utility for physically signifi- 
cant vibration problems, the case of the freely spinning 
shell with no hub has been treated only by the linear 
» method except for the boundary layer analysis. The major 
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result of this analysis is the demonstration of the continuous 
transition of the shell solutions to the flat plate solutions 
when the curvature approaches zero. 

By following the method of Reissner [13], Eqs. (5.2.1) 
and (5.2.2) can be reduced to a single differential equation 


* a Z 

/ c: 


(5.4.1) 


where 



nK(\-v*) 


X 


i 

>Ueb 


and 



Reissner obtains the solution to the homogeneous equation 
corresponding to Eq. (5.4.1) in the form 

■=. C f berjp + C x jj' + ^*5 Jf 
+C 4 kei if + Cg + e 7 

^ (5.4.2) 
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(5.4.3) 



Of the eight constants in Eqs. (5.4.2) and (5.4.3), Reissner 
shows that Cg can be excluded since only derivatives of F 
are significant and that must be excluded to preclude the 
possibility of circumferential displacement. 

A particular solution to Eq. (5.4.1) which has the 
proper limiting behavior as the radius of curvature becomes 
infinite can be obtained in the form 


u 9-4- Af - 


X a n r 


(5.4.4) 


By substituting Eq. (5.4.4) into Eq. (5.4.1), solving for 
the unknown coefficients, and equating real and imaginary 
parts, there results for the particular solutions 



W eo 


oo 


/ . ( -' r — 


o=o 



(5.4.5) 
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oO 




- l 


(5.4.6) 


The general solution to Eq. (5.4.1) is the sum of the 
homogeneous solution, given by Eqs. (5.4.2) and (5.4.3), 
and the particular solution: 


k»er £ + C? bet* j + C 3 kerj + C 4 kei j 4 C s 
•ajV jouJ-ft* V 7 (-/)* f£-f y ' 


(5.4.7) 




[ C, 6 e/ j - C^ker | +r 3 kei j - C+ ker£ 


'tizO-v'V 


8 ) 


The constants in these solutions are to be determined from 
the boundary conditions. We choose to exclude the possibility 
of a rigid body vertical displacement by setting C ^ = - Cf . 
The Kirchoff condition 
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V z t& - o s"t C - h 

leads directly to the conclusion that C 6 =0. The requirement 
of finite stresses and displacements at the origin is violated 
if and are not zero, since is singular at 

Hence Eqs. (5.4.7) and (5.4.8) simplify to 


= C, [ ber(j) - 1 ] +C Z >oei(j) 


4 





(5.4.9) 


. O-vW.4 ^fvWv+T 5 (-0”~L. (£t n 


(5.4.10) 


The two constants and C 2 can be determined from the 
remaining two boundary conditions 
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( b) 


tL O 


These two conditions lead to two simultaneous equations 
in the unknowns and C 2 * Before writing these two equations, 
note that at r=b we have 

L. _ k _ - L b * Et 1 * 

■f * J *r&. ~ 1 dz.*- J 

Y ^ 

From Eq. (2.2.1) with z=0 at r=b we find that 

-- 

By solving for R in terms of h and be we obtain 


R 


~z V, 


Within the scope of shallow shell theory this becomes 




•Z. Vi 


Hence we can write 
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J x tl'iCi-v' 7 ) /W - A 

(5.4.11) 


With this observation we write the two equations in and 


\ be / 1 \ 

jYrz(t^) *■ * A 




o 


7k[C,[t«f A -hW)^] -C t [ber A -(,-»>) ] 


“ U,"cA) + 7 - ’ o 


where 



■\+V pj** 1 r ^ 'y* Cdl H 
« „T 0 te^I z 



4n 



p-vO^P 


■a *2 


<oy 4 





and differentiations are denoted in the following way: 
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1/r - f'M 


a 

grg ker *4 * &er u } cAr. 


If we introduce the notation 


a„CX^ •* bei A 4 (/-vO 
a.-tCA”) - ber A - 0-v 1 ) 


feer'A 


bei *A 


«fe,M * ^ 

a,T.CX)- b -^- A 


^,a) - -$"(*)-•£ {/(A) 

- La) 


we obtain the set of equations 


s„c, - a,^c x ~ - J z <f>,(X) 


«*, c , - < 3 « c - 


/IRh^ 

B-i' 7 - 


4>, (A) 


which have the solution 


0/A) -- 


-g^) «W» - 
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J^\J 

/^(A) St, (A) - 4^60 <?i ( £A) 

a^CA) a, z ol) -3„(<1)3 ZZ (A) 


A slightly more compact form of the solution can be 
had if we note that 


„ , , x w p-j-jf i 4 . v' rf+'Kii] fA) 4 ” 

S 1 W) ' - ** *eo A L(— )ir U; 



We write these as 

-t«a) -- ^ J,M) 

s# -fco • -^V A.u> 


where 



2 


[(4.*l+0(4n+-s) + V^C4rl+^] 

L(*Z.>n+^) !j X 
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c-/rVr.*.) 

\£?.r>+t) Ip 



With these definitions the constants can be written 


C'C A) - <*6» 


where 


lp (\) - &n.CA) 

' " «hM)^A) - a M a)c^CA) 

<a tl CA) 

^CA') - <9 m (A) ^CA) 


Equations (5.4.9) and (5.4.10) together with the 


4.12) 


4.13) 


4.14) 


4.15) 
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constants defined by Eqs. (5.4.12), (5.4.13), (5.4.14), 

and (5.4.15), completely specify the stress and displacement 
distribution in the freely spinning shell. The stresses 
can be written explicitly by noting that since 




i jr 
7 jr 


M 


«© 


di r^ 


• - 7- 

- = ,o^V z 


f 


we can write, using Eqs. (5.4.9), (5.4.12), and (5.4.13), 



- ty(A) 


bei 'A(%) ,, x l *er 1 \(%) 

A (,<'/W\ " ^ ' A COt) 




| 


L-<r ‘(no) 

l T 



(5.4.16) 


'V * <P,(A)Cter A(^) - ~j- r)^ ] 


L<z. + x)!r utJ i <5 - 4 - i7) 

The normal displacement t<5” can be obtained by substituting 
Eqs. (5.4.12) and (5.4.13) in Eq. (5.4.9): 
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u3" “ 





«*) -')+ <VA) A(4)] 


+ 


J_ 


r xJL+rt c-ir 

\ W £-1 I ]* 



(5.4.18) 


Although not required for the dynamic problem, the 
calculation of the shell bending moments can be accomplished 
by noting that they are defined by 


M rr u)?* r + r* yr) 

^6© - - O ( y: r + V* c*^. r ) 


By performing the appropriate operations using Eq. (5.4.18) 
we obtain for the bending moments 

Vi(A ) 1 1*,- \c%) * o -vOb^v%)] 

-L iZ/£\4- (-/) * vv# +s ) * v^4w*fc)J/A 

UJ 4 IF 
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K/t 


QQ 


I? 



\oer a(oO 

/\c^) 


- v' l=>e/ A (%>)] 


+ \a ^ ^ H r /b)l 


(5.4.20) 


5.5 The Limiting Case of the Freely Spinning Flat Disk 
Consider the case when the radius of curvature 
becomes very large, or, alternatively when the parameter 
defined by 


X 


/ 


X - V ^4: 


where h/t is the ratio of the shell rise to the shell 
thickness, is very small. It can be verified that for 
small \ 


a,.W 

a.zCM 


Wei X + X z 

berX ^ 


9t.,( X) 


fee/ 'A 

~T~ 




-2 


- 114 - 



3 n a) - 


bee' A 


£ 

IC» 


i _4 y * 
4 


a z ,(AA 3,^X1 - a,/X) e xx C\) ~ 

' ^ [(•*.« +o!p v’-J 


:* - 


/ 9 Z 


A) ^ 

^>(A) a^CA') - ^ X (A) 3 n CX) ** tc" 

£,( A) a^(X') - s ./ A') ^er A 1 

It is seen then that 

h’^, •«. 

A-^-o 

xJIL x ° 

Furthermore/ for small A 

be/ A( / £) J_ 

H*) * * 


- 115 - 



ber ' X(D 

xr/b) 


- -U Aff 

bet \(1) * i(\ £f 

beC k(K>) - I 

By referring to Eqs. (5.4.16) and (5.4.17), it can 
thus be seen that 


/ /’ kvi 

A -*^o 





(5.5.1) 



I 4%^/ X\ T ? 

- -CT71 h) I 


(5.5.2) 


Equations (5.4.21) and (5.4.22) agree with the results of 
Timoshenko and Goodier [28] for the case of the freely 
spinning disk. Furthermore, by referring to Eqs. (5.4.18), 
(5.4.19), and (5.4.20) it can easily be verified that 
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I i wn 
A-+& 
Z -+00 


**• o 


(5.5.3) 


f i m 
A-»© 
f?-*o o 



O 


(5.5.4) 


/ / K* 1 

A“*© 

eo 



O 


(5.5.5) 


Hence, it is seen that the bending theory provides displace- 
ment and stress distributions which show a continuous 
transition to the flat disk results as and 

5.6 Computational Results for the Freely Spinning Shell 
The solutions for the direct stresses and normal 
displacement have been computed by direct numerical evalua- 
tion of Eqs. (5.4.16), (5.4.17) and (5.4.18). In addition, 
the outer surface bending stresses , defined by 


t <3Y g - 


t 


(5.6.1) 
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have been calculated by direct evaluation of Eqs. (5.4.19) 
and (5.4.20). It should be noted that Eqs. (5.6.1) and 

(5.6.2) are somewhat unconventional in that they have the 
dimensions of force per unit length as do the present 
conventional definitions of the direct stresses. This 
definition was adopted to allow a direct comparison of the 
bending stresses with the direct stresses. 

To furnish a basis for comparison, the corresponding 
results have been derived by using linear membrane theory. 
In the membrane case Eqs. (2.4.4) and (2.4.5) become 


.cL - QTi + vvnuDV 1 ' O 

(5.6.3) 




(5.6.4) 


The membrane solutions can readily be verified to be 


o 


N/ 


ft 


(5.6.5) 


(5.6.6) 


u> 


■VfV^ 

\ 


ra=^*H(Aa1 


(5.6.7) 


B rr 


<o 

VizO-v'*) 


0+\r) 

X x 


(•2, + V^)p 



(5.6.8) 


„ <> (' + *0 , v ■» , i 

6 ee" ~ 

(5.6.9) 

The results of the bending theory computations are 
shown in Figures 20 through 24. We have used several values 
of X varying from the nearly flat configuration to the 
shell with a fairly substantial curvature (or alternatively. 
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very thin) , The results shown were obtained with sf * o.to • 
However, Poisson's ratio does not contribute strongly to 
the character of the solutions. These results can readily 
be compared to the membrane theory of Eqs. (5.6.5) through 
(5.6.9) . 

The transition with decreasing A from a solution 
which is very close to the membrane results to one which 
practically coincides with the flat plate theory is clearly 
evident. It is noted that for high \ the bending stresses 
are small and show little variation except near the edge. 

By referring to Eqs. (5.6.8) and (5.6.9) it can be seen 
that membrane theory accurately predicts the magnitude of 
the bending stress in the interior of the disk but, of course, 
fails to do so at the edge since it is not accounted for 
in the boundary conditions . With decreasing A the bending 
stresses first increase, then decrease as the thickness, 
or bending, effect becomes predominant and then the decrease 
in initial curvature dominates . 

The direct stresses show a similar trend with decreasing 
A . For high A the direct radial stress is nearly zero 
while the tangential stress displays the parabolic character 
predicted by the membrane theory of Eq. (5.6.6). With 
decreasing X we see a monotonic transition of the membrane 
type solution to the flat plate solution. It is particularly 
interesting to note the radical redistribution of the 
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direct tangential stress as X is decreased. 

A similar observation regarding the deflection u>* is 
also made. The characteristic parabolic solution of the 
membrane solution is seen for high X, while zero normal 
deflection is approached for small X , in which case the 
inertia load tends only to stretch the midsurface. 

From these computations the importance of the bending 
contribution is seen and it is concluded that Johnson's 
results [8] may lead to serious errors. It is concluded that 
for values of X which are not large the full bending 
theory should be employed for computation of the equilibrium 
state even if a membrane theory is subsequently used for 
the vibration analysis. 

5.7 The Limiting Case of Large X -Existence of Boundary 

Layer Phenomenon at Outer Edge 

It is noted in Figures 20 to 24 that as X becomes 
large the direct stresses N and N QQ and the displacement 
approach those predicted from the membrane theory. The 
bending stresses and Bgg which are calculated from 

membrane theory are seen to approach those computed by the 
general theory for large X in the interior of the shell, 
but deviate significantly near the edge of the shell. In 
fact, we note that the membrane theory is incapable of 
satisfying the zero bending moment condition at the outer 
edge , and this accounts for the deviation between the two 
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solutions in this region. Since for large X the deviation 
between the general theory and the membrane theory is 
localized near the outer edge, a boundary layer phenomenon 
is suggested. Reissner [14, 38] has dealt extensively with 
problems of this type for laterally loaded shallow spherical 
shells. It is the purpose of this section to demonstrate the 
existence of this phenomenon for the case of the spinning 
shell. 

Because of the ease with which a nonlinear membrane 
solution can be obtained, it is convenient to investigate 
this phenomenon by using Reissner' s nonlinear equations for 
the spinning shell given by Eqs. (2.4.8) and (2.4.9) and 
repeated here 


D 



+ 


7 

r ~ r*- j 



o 


(5.7.1) 





(5.7.2) 


where 


r*lrr ' P 

(5.7.3) 
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The boundary conditions are that the solution be regular 


and symmetric at the origin and that at the outer edge, 
r=b 



It is convenient to non-dimensionalize equations 
(5.7.1) and (5.7.2) by introducing 

S' $ 1 

4 > % 4>o j (x) 

in which case we obtain 





* O 


p t ^ - £]- £*? **(*)£■£ -'gppJ&h 


If we let 


v* /^/y- &* 

Et'bJ "' <? 


*.«£X*XW>“ le *) ' (>Xi) v' 




where 


if -- 


W 1 i 9 ^ 


then we have 




( 5 . 7 . 6 ) 
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( 5 . 7 . 7 ) 




•«- - X 


It is seen that ^ is a measure of the nonlinearity of the 
problem. 

The stresses and displacements will become 




f f 




Q& 


~"io $'(*) + '^( r ) 


ay 





By expanding and using shallow shell relations and the 
definition of X we obtain for 
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duS. J_ 

£ V*To^v*) 



The bending moments will then become 



The boundary conditions on the non-dimensional form of the 
problem will be that f(x) and g(x) be regular and symmetric 
at x=0 and that at x=l 


1 




o 


six') "*■ o 


Note in Eq. (5.7.6) that 

1 && ^ > Jt_ J_ 

lo 4 et ^ ^ " \ 4 


Hence, we write the defining equations 
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~ o> 


(5.7.8) 


*x* + X - + 4f 

(5.7.9) 

To investigate the boundary layer phenomenon which 
exists for large X we make use of some of the concepts of 
perturbation theory [39, 40]. Since the outer edge is the 
region of primary interest we will consider "inner" and 
"outer" expansions relative to that point. A suitable outer 
variable is X with I t in which case Eqs. (5.7.8) and 

(5.7.9) become 

^l^xx* X + * ' -• ° 

(5.7.10) 

**x + 7 *x -*>■$ + 4 f ' 

(5.7.11) 

where is assumed to be a small parameter defined by 
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The first term in ain asymptotic series in powexs of £ will 
be defined by the equations 


* - 



O 


( 5 . 7 . 12 ) 


$ 


xy 





( 5 . 7 . 13 ) 


A solution to this degenerate set of equations which 
satisfies the boundary conditions at x=0 is taken as 


( 5 . 7 . 14 ) 

' U'-^l 

( 5 . 7 . 15 ) 

This is the nonlinear membrane solution in the case when 
the membrane shell is not deflected to a flat plate. We 
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will assume in this analysis that "Z V ^ I , the condition that 
the deformed membrane shell is not flat, so that the above 
form of outer solution is valid. It is noted that Eqs. (5.7.14) 
and (5.7.15) will not satisfy the boundary condition of 
zero moment at x=l, since it is found that 

.wU ' - eC'WK 

where 

& x \ - 

To investigate the nature of the solution in the 
neighborhood of the free edge at x=l it proves convenient 
to recast the problem by defining 

Sex') -• 

In this case Eqs. (5.7.8) and (5.7.9) become 

l 4 *X ~ + C * “ ° 

V ' (5.7.16) 

(5.7.17) 
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The new boundary conditions require that F and G be 
regular and symmetric at x=0, and that at the free edge 



- 0 + 




F" o 

The inner problem is formulated by choosing the 
parameter to be 


in which case 


X - 


X. <L 

£ b 


The outer edge is made the origin by the change of variable 

5 - j[ - -^-(v - V*) 

Equations (5.7.16) and (5.7.17) then become 



JL 









'-*1 


o 

(5.7.18) 






:nzF - C i - ^7 + 


V. ^ 

z i-£f “ ° 

(5.7.19) 
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The boundary conditions are 


1 


1 


o : 




F ' o 


OW) 


■s' 


and G regular and symmetric 


If ^ is 0(1) , that is if we restrict ourselves to a 
narrow region in which \ - c /\o is 0(6), then Eqs. (5.7.18) 
and (5.7.19) and the boundary conditions can be approximated 
to within terms of Ofe) by 
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(< -S)p- 


o 


(5.7.20) 


1 - 
1- 


d Z F _ 

- O 

o \ 

- 

C ' + m0 ^ 
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- oo t 
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(5.7.21) 
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Equations (5.7.20) and (5.7.21) can be construed as 
defining the term of 0 (€°) in an asymptotic series in powers 
of the small parameter € . In principle, the higher order 
terms could be sequentially obtained and the procedure of 
matching with the outer solution could be accomplished. 
However, we have succeeded in demonstrating that there does 
in fact exist a narrow region in which the solution departs 
radically from the outer solution and can be made to 
satisfy the free edge condition. 

From the above analysis and the behavior shown by the 
computed solutions, it is concluded that for large A the 
solution is characteristic of membrane theory over the 
interior of the shell but in a narrow region at the edge it 
will involve significant bending to allow satisfaction of 
the free edge condition. 

A complete analysis of the boundary layer phenomenon 
will be undertaken as an extension of the current effort. 

5 . 8 Linear Theory for the Case of a Fully Clamped Central 

Hub 

In the case of a fully clamped central hub the boundary 
conditions to be satisfied by the general solution given by 
Eqs. (5.4.7) and (5.4.8) are 
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(5.8.1) 


(5.8.2) 


The Kirchoff condition again requires that Cg vanish. The 
general solution then becomes 
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C, J 4 


Cj be/ J; t /rer j jf -f + C«s 


+ +V*) C - 1 ) 


Z 1 C-0 y ’ /rV 


4 -*l 4 (© 


(5.8.3) 


F * - \^C» bet ^ - C-^ ber^ + C-^ ktf/ - C4 ^ 


i-V/Xs4 "tOM ^ (-o’"" /-£\ 

^ U; + x* sL, ^ 


va " 1 - r \ 4 vi + A- 

(5.8.4) 


Me* -- (77^) [ C i bei ’(£) - C-l <=>«■* '(X) + C s te.' '( 4 ) 

-C<kcr'($) J ♦^f.C*) 


(5.8.5) 
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The constants C^ through C,. are determined from the boundary- 
conditions which lead to the following five algebraic 
equations (here written in the same order as the boundary 
conditions of Eqs. (5.8.1) and (5.8.2)): 

■*• - C, her \(%) - C z hei ~ C-^ker \(^/io) 

-c* kc; A (%) - ' zJ ^ L) ^(4 (5 . 8 . 10) 


+ke /A^)] 


x («/o 


-V C^r x< 3 4) - + + ke,-\(^)] 


(5.8.11) 
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(5.8.12) 


C, t>ei' A -<\be/A +C i U«i'A - C 4 ker'h «• - 

(5.8.13) 

- C*\ ber A - ] 

+ C^,\wei A * <"^] - C^ker A - <■'-*> ^ 1 

_ "X + V'* q ( (5.8.14) 

^T. (Kl^/ 

It is seen that through are defined by the four 
Eqs. (5.8.11) through (5.8.14) while is obtained from 
Eq. (5.8.10). As soon as the constants are determined, the 
direct stresses, bending moments, and normal deflection can 
be computed from Eqs. (5.8.5), (5.8.6), (5.8.7), (5.8.8), 

and (5.8.3), respectively. The outer surface bending stresses, 
defined by Eqs. (5.6.1) and (5.6.2) are directly obtainable 
from the bending stress results. 

For demonstration purposes we have evaluated the stresses 
and displacements for the case of a hub with ( a /b) ■*- o.t-j-s 
and XT’ - o.-2o . The results, plotted for several values 

of \ , are given in Figures 25 through 40. The character^ 
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istics of the solution will be discussed in Section (5.10) 
where they are compared with the nonlinear results. 

5.9 Nonlinear Theory for the Case of a Fully Clamped Central 

Hub 

As will be seen in detail in Chapter 6 , the nature of 
the spinning shell vibration problem is critically 
dependent on the stress and displacement distributions 
obtained from the equilibrium solutions. To assure that 
an accurate evaluation of the equilibrium state is available 
for the dynamic problem, extensive calculations have been 
performed, for the case of the clamped central hub, using 
Reissner's nonlinear theory [14] which allows finite 
rotations of shell elements. The linear and nonlinear results 
have been compared and the regions of parameter values in 
which the nonlinear effects are important have been established. 

The governing equations of Reissner's theory were 
derived from a variational principle in Chapter 2, were used 
in the boundary layer analysis in Section (5.7), and are 
repeated here for reference: 

=<=> 

(5.9.1) 
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r ti fr -- i/y 


(5.9.3) 



c/ (4/ 

c Jr 



LO 


V z 


(5.9.4) 
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(5.9.5) 


As shown in Section 5.7 , a non-dimensional form of 


these equations can be written 


■f + is. - JL-= ♦ ^ v -^ 1 • - X 

** x x* ff ^ v 
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The functions f and g are defined such that 


M fr =. (T,4v')pui‘\o’ 1 ~- 

M^-C^rtpufb'lVW 

_ c^+v*) „ z„7. - N 

*** " Et f? P^ ^ 


(5.9.8) 


(5.9.9) 


(5.9.10) 


- 140 - 



Equations (5.9.8), (5.9.9), and (5.9.10) can be written 

in the slightly modified form 


- 


_ C^4vQptO t l? Z g/.b.V 2 

t " £ £ V X 




(5.9.11) 






= Ve 


(5.9.12) 


K 4jy> 

Vo 2 - ci y 




- tf^x) 

where for convenience we have defined 


V>(y) -- $'<.*) + 


jr 

■% + V' 


(5.9.13) 


The definitions of Eqs . (5.9.11), (5.9.12), and (5.9.13) 

prove to be convenient when describing the equilibrium direct 
stress and displacement distributions for the vibration 
problem in Chapter 6. 

The bending moments are defined by 
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(5.9.14) 
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(5.9.15) 


The outermost surface bending stresses are given by (13] 


- 

(o 

ft 

**■ 


+ k ^©© 
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By referring to Eqs. (5.9.14) and (5.9.15), and making use 
of the shallow shell approximations, it is found that 




i fizU^FV ' L «* 


+ \T 


f] 


(5.9.16) 








(5.9.17) 


The boundary conditions accompanying Eqs. (5.9.6) and 
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r 


(5.9.7) are 



(5.9.18) 


(5.9.19) 


Equations (5.9.6), (5.9.7), (5.9.18) and (5.9.19) 

constitute a nonlinear two-point boundary value problem. 

A method which is essentially a modification of Archer's 
technique [18] has been employed to obtain solutions for the 
stresses and axial component of displacement. 

Archer's method is iterative in nature. He solves the 
linear portion of the problem and then utilizes this solution 
to evaluate the nonlinear terms. These are then treated 
as an additional nonhomogeneous contribution to the linear 
equations for the next solution sfep. This procedure is 
carried out for several iterations until successive solutions 
coincide to within a specified error. Archer's chief 
contribution is the finite difference scheme which he has 
found to be appropriate for solutions to the successive 
linear two point boundary value problems. 
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Archer's method, in its basic form, was initially 
employed to obtain solutions to the present problem. The 
grid employed 36 points in the interval from to 

X l and represented nearly the limit in complexity 
which could be handled on the Wichita State University 
IBM 1620 digital computer. While some results obtained by 
this method were satisfactory, several instances occurred 
in which the iteration procedure diverged. Since a buckling 
phenomenon is not expected, it was concluded that the lack 
of convergence was due to some inadequacy in the technique. 
No attempt was made to identify the source of the conver- 
gence problem and it is possible that the particular 
implementation of the method was not adequate. 

In view of the convergence problems encountered with 
Archer's basic technique, a slightly modified method was 
employed and convergence was obtained for all combinations 
of parameters used. This method, discussed more completely 
in Reference [41] , employs the finite difference method of 
Archer, but defines the successive linear two point 
boundary value problems involved in the iteration in a 
different way. The method is best explained by considering 
a simple example. Assume that we have the simultaneous 
nonlinear equations 
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(5.9.20) 


4" + m -- a(x) + 4 ? 
ip + 4> x 4* 


The first step in the procedure would be to obtain 4i , 
and 0, , solutions to the linear problem 

4,'* + 4, -- a ( y) 

l/> ( + <4, ' ° (5.9.21) 

We then write the true solution as the sum of 
and some remainder terms and 

4 > 4 , + 

0 ^Pt f (5.9.22) 

By substituting Eqs. (5.9.22) into Eqs . (5.9.20) and making 

use of Eqs. (5.9.21) we obtain the governing equations for 
4 t and : 

4^" + 4-z. - 4 , ■- 4, z + 4^ 

+ 4^ - - 0, 4^ x 0, 0, + 0> z 

(5.9.23) 
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In Eqs. (5.9.23) , 0, and are known functions of 

x . 

We now repeat the above procedure by considering 0>-^ f 
l/j ^ to be the solution to the linear portion of Eq. (5.9.23) 

0* + ^ - 0, z 


0 + 0s - 0, A - <A 0* ' 0, 0/ 


(5.9.24) 


We then consider <k' A to be composed of and 

some remainder 4*. *4 ' 

<t>\ ' <£ % + ^4. 

' 0V + 4>4 

^ * * (5.9.26) 

The equations for <t>4 . 0* are obtained by substituting 

Eqs. (5.9.26) in Eqs. (5.9.23) and employing Eqs. (5.9.24): 


<£>4 + 04 - “*■ 0 , 04 ” 0 * 04 ~ 0 * + 4>4 

04 ,f + 04 ' 0, 04 “ 0, 04 “ 0 3 04 - 03 04 ' 03 03 + 04 0/ 

(5.9.27) 

In these equations - ^3 are known functions . 

We proceed as before, solving the linear part of Eqs. (5.9.27) 
and proceed through n such steps until I I and 
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become small compared to the accumulated 


n 

4) oz 4 > ( + ~ Z, 

i “ / 
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tp 'a ip t + ip z + • • ~ Z, c// xi _, 

t e / 

The main advantage of this method is that at each 
iteration after the first some features of the nonlinearities 
are included in the linear equations, while in Archer's 
method the particular form of the nonlinearity is relatively 
unimportant at a given solution step. The principal dis- 
advantage is, of course, the increased complexity and logic 
involved. It can be seen from Eqs. (5.9.21), (5.9.24) and 

(5.9.27) that the linear equations to be solved change 
from step to step. However, the change is systematic and 
can be accounted for simply by updating the coefficient 
matrices in the finite difference scheme after an iteration. 
Other than the necessity of a few more computer instructions, 
the computational requirements for this method are not 
substantially more severe than in Archer's basic technique. 

In the application of the above technique to Eqs. 
(5.9.6), (5.9.7), (5.9.18), and (5.9.19) convergence was 

obtained for all the combinations of parameters which were 
tried. The results of these computations are plotted on 
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Figures 25 through 43 where they are compared with the 

3 . 

corresponding linear results for g- = 0.125. 

5.10 Comparison of Results for the Linear and Nonlinear 
Theories 

Solutions for the stresses and displacements in the 
spinning shell fully clamped at the hub have been obtained 
by using the linear theory of Section 5.8 and the non- 
linear theory of Section 5.9 with a grid of 36 points in 
the region o.rz-s< X ^ I .We were required to impose 
the restriction of 36 grid points because of computational 
limitations. In some of the results some slight deviations 
from more exact results will be noted but the general 
agreement between comparable calculations is considered 
satisfactory. The results have been plotted in Figures 
25 through 43 for four values of the shell geometry 
parameter X - o . o-s*&> , /.o , ■& . o , and “7.0 , 

which cover the case of a nearly flat shell to one with 
a fairly substantial curvature. For each value of X we 
have obtained the linear solution and the nonlinear solution 
for five values of ^ , which is the inertia loading 

parameter and scales the influence of the nonlinear terns. 

We have taken V -0.0 1 , , /.O , ^.o , and 

12.0. O f which cover the range from very low to very high 
inertia loading. A value of "s-o.S would correspond to 
deflection of the shell to a flat plate in the nonlinear 
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membrane shell theory (see Eg. (5.7.15)). 

Figures 25 through 28 show the radial direct stress. 

For small A it is seen that the stress distribution is 
nearly that of the spinning flat fully clamped disk. 
Furthermore, for X *. o.osftO and practically no 

nonlinear effect occurs, as is seen by the fact that no 
variation with V occurs , to within the accuracy of the 
graphs, and by the fact that the linear and nonlinear 
theories are in good agreement. For increasing A we 
find a substantial variation of the solution with Y • For 
low Y the nonlinear and linear solutions nearly coincide 
while for large Y the solution is nearly that of the flat 
disk. In particular, for small Y and increasing X , 
the characteristic decrease in importance of the radial 
stress is noted. 

These phenomena are readily explainable on physical 
grounds. For small X the shell is very nearly flat so that 
we expect the stress distribution to be near that of the flat 
disk. For higher X curvature effects become important 
and for small Y a membrane stress distribution is approached. 
For increasing Y the disk progressively flattens out and 
the stress distribution alters radically until finally, for 
very high X , the stress distribution is about the same as 
if the shell had been flat to start with. In this case the 
stress distribution required to overcome the initial curva- 



ture is small compared to the additional stress built up 
after the disk is essentially flat. 

Figures 29 through 32 present computed results for the 
tangential direct stress. The same trends are seen here as 
were seen in the case of the radial direct stress. For 
small \ the stresses are those which occur in a flat 
disk and no nonlinear effect is possible. For high \ and 
small the linear membrane results are approached, and 
for increasing ^ the nonlinear effects cause a radical 
variation in the stress. For very high the flat disk 
results again are approached. 

Figures 33 through 40 are plots of the magnitude of 
the maximum bending stresses. The bending stresses were 
computed from the relations 
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The non-dimensional forms plotted in the figures, are 
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For X'-O.O'S^j , the shell is nearly flat, and practically 
no bending stress is generated, and no nonlinear effect is 
seen. For A-/-0 r the shell is still relatively flat so 
that very little nonlinear effect is seen until ^ becomes 
very large and the disk is effectively flattened out. The 
bending stress shown for ^ -1*20 would be essentially the 
residual bending stress required to flatten the shell and 
would show almost no variation for further increases in 
^ . The parameter ^ could also be lowered considerably 
before significant variation in the bending stress is 
seen. When X is increased to 3 and above the nonlinear 
effects appear for relatively small values of ^ . The 

effect of increasing is to decrease the bending stresses. 
For high V the residual bending stress becomes small 
for the higher values of X as the membrane effects become 
predominant. The characteristic initial increase in 
bending stress with increasing X up to a point, and then 
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decreasing bending stress with increasing X is clearly 
evident . 

Figures 41 through 43 present the results of the 
nonlinear theory for the bending slope. We have not 
included the case As 0. 05<£since the results for the bending 

_7 

slope are on the order of 10 and the single precision 

accuracy is open to question. For the case O we 

have used a scale factor of 10 in order to maintain a 
uniformity in the figures. Further, note that we have 
employed the nondimensionalization 



which can be written 


c# 



z±v ^ /— ) 


A particularly interesting feature occurs for large X 
when If is large. Note when $*IZO with X r 3. 0 or 
7 . 0 the slope of the deflection curve is nearly 
linear. In this situation the disk has practically 
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flattened out and the bending slope has become the negative 
of the original shell slope. This effect is not nearly 
as predominant for smaller X . This implies Vt is 
small and both h small and t large would reduce the tendency 
to flatten out. 

The result of primary interest from these computations 
is the strong nonlinear behavior of the solutions for 
increasing ^ r particularly when X is not very small. 

The substantial variation of the solutions with changing 
suggests that considerable care should be exercised if 
anything but the nonlinear theory is to be employed for the 
computation of the equilibrium configuration. 

5 . 11 Conclusions Regarding the Appropriate Theory for 
Computation of the Equilibrium State 

As will be seen in Chapter 6 , the equations of motion 
of the spinning shallow spherical shell for free vibration 
about the equilibrium configuration depend strongly on 
the direct stresses and displacements which exist in the 
equilibrium state. For this reason particular care should 
be exercised in choosing the appropriate theory for the 
analysis of the steadily spinning shell. 

By referring to Figures 20 through 43 and the preceed- 
ing analyses it is seen that for large X a membrane theory 
is very nearly exact except for narrow regions near the 
boundaries where bending is important. The width of these 
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regions decreases as X increases* 

For smaller values of \ , bending becomes important 

and a continuous transition to the flat plate results 
( requires that it be considered. 

Small values of the inertia loading parameter 
give rise to negligible nonlinear effects but these effects 
become increasingly more important as ^ is increased. 

The effect of increasing is more predominant when X is 
large. 

It is concluded that any general purpose vibration 
analysis should be based on the nonlinear theory which 
includes bending effects. The use of a membrane theory, 
as employed by Johnson [8] , is not appropriate, particularly 
for small X • In addition, the effect of bending in the 
edge zones on the vibration characteristics should be 
examined in the membrane case. The use of a linear theory 
may well be adequate for most practical applications but 
its use would not substantially simplify the vibration 
analysis . 
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CHAPTER 6 


THE EQUATIONS OF MOTION FOR THE FREE 
VIBRATIONS OF A SPINNING SHALLOW SPHERICAL SHELL 

6 . 1 Introduction 

The study of the transverse vibrations of a spinning 
shallow spherical shell presents two interesting features 
not found in analyses of non-spinning shells. The first, 
and most significant, variation introduced by the presence 
of spin is the generation of a pre-stressed equilibrium 
configuration about which the vibrations occur. This 
configuration has been studied in detail in Chapter 5. The 
second new feature which arises is the presence of a 
coriolis coupling between the three deflection components. 
These two additional considerations combine to produce 
a mathematical problem which is different and considerably 
more complex than the corresponding problem for the 
stationary shell. It is the purpose of this chapter to 
fully formulate the problem for numerical computation by 
an established technique. The available computation 
facility has proved to be inadequate to handle a problem of 
this size so that no numerical results have been obtained. 
Efforts are under way to pursue the matter further at another 
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facility. 

6.2 Basic Equations — Meridional, Tangential, and Normal 

Deflections 

The differential equations which goverii the small 
vibrations of a spinning shallow spherical shell about its 
equilibrium configuration were given in Eqs . (2.5.9) , 

(2.5.10), and (2.5.11) for the case when the deflections 
are resolved in the meridional, tangential, and normal 
directions relative to the undeformed shell. They are 
repeated here for convenience with the now unstarred 
quantities understood to signify perturbation quantities 
and the quantities subscripted with ”0" understood to 
represent equilibrium values: 


( 6 . 2 . 1 ) 
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( 6 . 2 . 2 ) 
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v T v 7 u 3 - - 7T^ (^07. tS-a r + r^> o u> r ) 

" i>i> (" r ^© 0 * u ^v + ^e) 

= y 03- --Z .UJO-^ -(u4 )cO^] 

(6.2.3) 

In addition to the equations for dynamic equilibrium we 
must include the mid-surface strain-displacement relations 
for the perturbation quantities, here assumed smalli 
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+ u>o r U> r 


(6.2.4) 
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(6.2.5) 
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( 6 . 2 . 6 ) 

To complete the specification of the problem we state the 
boundary conditions for a shell of outer radius b with a 
fully clamped hub of radius a: 
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(6.2.7) 

r -- t» : <gr r •«- o 

° 

•- -d(^ + 7«>,) -- o 

hCi&CHorfr CoND: ^:V Z uX+ ~r~ ( T 7 = 0 (6.2.8) 

For the free vibration problem at hand we make use of 
the results of investigations of stationary shells to assume 
the following solution forms: 

uCr,©^)- U(r) Cos f 

i>(r,e ; 4) - Vcr) hg + tO-0 

* V1/^ CoS(/1©+77Vj (6.2.9) 

With this form for the displacements the strains are given 

by 
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4 r =(u r + ^ * “% r W r ) cosfn 0 + -pi) 

^-(■F^-F^) *<>*<"* **0 

V fs * (V r - ^ - '^r - ” uJ °r'2r)^n(r'e+i3i) 

( 6 . 2 . 10 ) 

By utilizing the stress strain relations 

'S'r --F^[ e r + v'£ s ] 

= 7 TF>-l 6 e * v ' £ f] 

r ra * G, X fe ( 6 . 2 . 11 ) 

we find that the stresses can be written 
<T r - S r Co’S ( /)© + *? ^ ) 

<^r e - S e co-^c *© + ~p-l) 

= S re /l e + ) (6.2.12) 

where 
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The equations of motion can be written 


(6.2.13) 


^rOS,.) - S e + ^S ro = -^^UfZcof)V + (u + |- r )^] 


(6.2.14) 


57^5 re )- M S 6 tS rc - -yy<r[-fi l V+-Zu^-pCu + '^') + ^ z v] 

(6.2.15) 

f L 4 w/ -?i(ru> 0r S r + rcr^ o <) 

4 r (7 ^ c \A/ - vn o>o r ^ce) + $ e ) 

- mV^W + 'Z«D 1 0\/-£ +(U<- W £) Ui Z -~ ] (6.2.16) 

where 

,^.il id. 

L ' dr 1 r dr ~ ~pi 

By using Eqs . (6.2.13), Eqs . (6.2.14), (6.2.15), and (6.2.16) 

can be written 
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It proves convenient to cast Eqs. (6.2.17), (6.2.18) 
and (6.2.19) in non-dimensional form by introducing 

r 
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U u £v) 
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Further, from Chapter 5, Eqs. (5.9.11), (5.9.12), and 
(5.9.13) we have in the equilibrium configuration 

— ^ ■«. i £ ujZ !? Z 
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(6.2.27) 


By introducing Eqs. (6.2.20) into Eqs. (6.2.13) we find 
that the stresses in terms of the displacements can be 
written 


o _ i-v" £ WJA _U_^-vW£ u ^ 0 N )~r( > 

'Ve - "2 f-\r* \ ^vr X X '•fe?' y/ y j 

With Eq. (6.2.27) we can then obtain 
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Equations (6.2.29) can then be written in the convenient 
form 
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where / S re are defined by the bracketed 

terms in Eqs. (6.2.29): 

S r “ Uy + u3- + a> x V- + t*> *- vr ^) 

Sq — ^ +• u3- + *■£ + cd-y) 

% ce ~ °l( - if “ ~ ^ y" (6.2.30) 

By utilizing the relations of Eqs. (6.2.20) through (6.2.30) 
we obtain the non-dimensional equations of motion 

y u xx + * + u x + 0- t +\r)x ■+0-V') 

+ V* fc(r)R -1 . +*Vv)]# - 

+ - O (6. 2 .3i) 

^ - *u&« t f ^ - ^V) «>« - ^ f 

+ 44? ^ (^) Z< ^\.' + (^) T ]*> + "*(wXu+<u>)[ - O (6.2.32) 

Cu>- +f(x)<Ar] 

+ * -ifhvl ?wv»] + 

(6 - 2 ' 33) 


- 164 - 


The corresponding non'-'dimensional boundary conditions are 
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It is seen that the non-dimensional equations of motion 


depend on four parameters: 
(■%+y’) w\ ui 1 1? 1 
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In addition, the boundary conditions require the specifica- 
tion of the ratio of the hub radius to the disk radius. 

It is observed that the dynamic problem requires the 
specification of one additional parameter as compared to 
the equilibrium problem. It can be seen by referring to 
Eqs. (6.2.31), (6.2.32), and (6.2.33) that the extra 

parameter, ^ , scales the relative importance of the inertia 
terms in Eqs. (6.2.31) and (6.2.32) and the Coriolis and 
centripetal acceleration terms in Eq. (6.2.33). Since the 
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shallow shell theory used here implies the neglect of terms 
of the order of it i s concluded that terms of 

this order in the equations of motion can also be dropped, 
provided that only the transverse motions are of interest. 
This is in agreement with the results of Reissner [22] 
regarding the neglect of longitudinal inertia for primarily 
transverse vibrations of shallow shells and extends the 
conclusion to include the coriolis and centripetal coupling 
terms in the equation for transverse motion. As shown by 
Reissner, if advantage is taken of the neglect of these 
higher order terms the governing equations can be 
simplified by the introduction of a stress function which 
identically satisfies Eqs. (6.2.1) and (6.2.2) in the 
case when the right hand sides are zero. The equation 
satisfied by the stress function is found by employing a 
suitable compatibility equation. 

6.3 Basic Equations — Radial, Tangential, and Axial 

Deflections 

For some purposes it may prove more convenient to have 
the governing equations resolved in the radial, tangential, 
and axial directions as was done by Johnson [8] in the mem- 
brane case. The equations for this case were originally 
given by Eqs. (2.6.15), (2.6.16), and (2.6.19) and are 

repeated here with the stars and bars suppressed 
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( 6 . 3 . 1 ) 


^ ^re) +• + 7^ e - tTuJ w - o-co 1 

( 6 . 3 . 2 ) 

't v ^ ZuV - *• «v)^r + r<3 7 C) £ A'( > 

+ «**,.) ^ * 7 ^ o <*> e | ( 6 . 3 . 3 ) 

The corresponding strain-displacement relations for the small 
perturbation quantities are 


6 r ^ u r + u> v «%. - £ 

( 6 . 3 . 4 ) 





( 6 . 3 . 5 ) 



+ o- - j. «^ r c<3© 
f r r — 7 - - ^ 


(6.3.6) 
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The boundary conditions for the case of a fully clamped hub 
of radius a and a free outer edge will be given by Eqs. 

(6.2.7) and (6.2.8). 

By utilizing the wave type solutions given by Eqs. 
(6.2.9) the strain displacement relations can be written 


(6.3.7) 

- [ 7 T + n -pr ) Co ^ C * & + ■P’t) 

(6.3.8) 

"( V r ~ "7 - T * na \-7 

(6.3.9) 


The stress-strain relations can then be employed to yield 
^ C.o<s G 

(\s\& 

(6.3.10) 


"*■ ^ ^iki ^ ►'I Q + '^Q'i ^ 


where 



( 6 . 3 . IX) 


V - u, * ^ r w r -£w r . *($ * <$)] 

S e -TTr-l'T + T- + S ( U r * u>* r W r - ■£ W r )J 

In terms of the assumed form of solution the equations 
of motion can be written 

- - - ^r[(-p Z +uJ z ) U + ~ZoOf)V] 

(6.3.12) 

dr ( f ^r&) ~ ^ + ^>r© ~ -mr^(-f) l + uj^V + ~Zuj-p L/J 

(6.3.13) 

f^-7^W‘£ + “VK ^V'v'rf 

f(-^ + u> 0f )S re - (6 ^ 14) 

where 

1 7 . __ JL d_ yl 

L d^*- r dr ' 

The equations of motion in terms of displacements become 
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(6.3.15) 


r U fr , + f( u>c r - r 4?) W^ r + u r -t J. j(, -vXu>i f - r /g) t r(u> 0 - /g)~j V/ r 

- 1 1 + ^ c^e>] u - 

+ ^- - v ^ )vv1 ^ V [[,-*^/? U + - o 

(.-v-)rV' r - ^U r 

+ (izjrV^rjit + tgflv 4 i(£)uj - o 


(6.3.16) 


X i # M/ - r Jr*f r ("? + “V) r ®? 0 j* 

f r l ' r «o Vv/ - fuJi r) 3r# ' "« (6 3 _ 

where and S'fg are defined by Eqs. (6.3.11). We 

obtain non-dimensional equations in the same way as in 
Section 6.2 by employing Eqs. (6.2.20) through (6.2.27) 
together with the appropriate non-dimensional form of the 
stress-displacement relations for this case: 

S,- ' JTTi u x + [ ~ y ] Wjr + ^(7 + 

$ e |t + TT + Wtf* +1^6) -x] u>J[ 

S, e - ^ ^ - f + "[*- *07 f 

By introducing the non-dimensional stresses S r r , 

^C© r defined by 
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V ' + [ VjW - *] + v'C % + 

- ‘y + \T ( Uy + ( ^6() “*J ^x) 

»CB' ** ' +4*- ^J6f)] -f 


we have for the dimensional stresses 
o E \a/q 

'V -- ,. r z, ^ 

s, 
s. 


£ Wo 
>© ' /-v'' z - R’ 

i-r e_ w? - 

y© ' /-r*- R <e 

The non-dimensional equations of motion for this case are 

Y U xy + *t #%(*)“*]“%* + ^sgp^ O-rfoj/ffl-xJ+rffyxJ-xkfox 

-{, * „* u^>] a. - *f« -*J f 

W*W - o 


(6.3.19) 


LlX'vo- *(t +r) 
** * *y ~ ■ -z _ 


u * „ x ] u> x 


+(*y*) -*)]-* •‘■^ o 

1 7 (6.3.20) 

l }«> + ***)«>*[ 


_ A 4 * 1 /d£p 


r© -s+v^ 


(6.3.21) 
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Equations (6. 2 . 34) and (6.2.35) can also be employed to 
specify the boundary conditions for this case except that 
the appropriate definitions of and must be 

employed in Eqs. (6.2.35). As in the case of meridional, 
tangential, and normal deflections we find that the 
differential equations and boundary conditions depend on 
five parameters. In addition, we once again note the 
relative unimportance of the longitudinal inertia terms 
in the membrane equilibrium equations. 

There does not appear to be any particular advantage in 
choosing one or the other of the two forms of the governing 
equations of motion for purposes of numerical computations. 
Although the form valid for radial, tangential, and axial 
deflections appears to be slightly simpler in the bending 
equation, a few manipulations can be employed to reduce the 
other equations to a nearly equivalent form. 

6 . 4 Numerical Computation of the Natural Frequencies and 

Mode Shapes 

The method which is proposed here for the numerical 
computation of the natural frequencies and mode shapes is 
an adaptation of a technique reported by Zarghamee and 
Robinson [26] . Reduced to its basic concepts this method 
is an extension of the technique used in Chapters 4 and 5 
to determine the natural frequencies for flat membrane 
disks. 
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Since the problem is linear we know that a solution 
can be obtained as the sum of a number of independent 
fundamental solutions. In particular, we choose to define 
our fundamental solutions as those solutions which satisfy 
the boundary conditions at one boundary but which have 
unspecified values at the other boundary. By superposition 
of these solutions with suitable arbitrary constants we 
obtain a general' solution whose constants can be determined 
by the requirement that the remaining boundary conditions 
be satisfied. This procedure leads to a set of homogeneous 
algebraic equations for the undetermined constants which 
can only be satisfied in a non-trivial way if the 
determinant of the coefficients vanishes. The determinant 
of the coefficients is formed from the fundamental solutions 
and certain of their derivatives evaluated at the boundary 
in question. Since these solutions depend on the undeter- 
mined natural frequency we find that the problem reduces to 
that of finding values of the natural frequency for which 
the determinant of coefficients vanishes. 

By referring to Chapters 4 and 5, which discuss the 
vibrations of spinning membrane disks, it can be verified 
that the fundamental solution was taken as that solution 
which satisfies the finiteness condition at the outer edge 
of the disk. The vanishing of the deflection at the hub 
leads to the condition that the fundamental solution itself 
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vanish at the hub. It is seen that the determinant of 
coefficients in this case reduces to a single element. 

The disk vibration analyses led to exceptionally simple 
eigenvalue problems since the boundary conditions were 
simple and the form of the fundamental solutions was known 
apriori.. 

The problem at hand is conceptually the same but 
practically considerably more difficult. The boundary 
conditions, of course, are more involved, requiring four 
conditions at each boundary. However, the biggest complica- 
tion is that the form of the fundamental solutions is not 
known apriori, but rather, must be determined by numerical 
integration of the equations of motion. The eigenvalue 
equation requires the determination of the zeroes of the 
determinant of coefficients whose elements are obtained 
by numerically integrating the equations of motion. 

In the present problem we will identify four linearly 
independent fundamental solutions which satisfy initial 
value problems defined as follows s 
Cl) 0(u, ~ o 


«(%) 

U k (*4>) - ' 


u-(aA) ~-o 



u>6=*/b) * ° 

u± KX (a/b) ~o 


u%(a/b) * o 

<Axx(*0 ’■<=» 

(6.4.1) 


Yields fundamental solution U,(x) t (x) , U?, (x) 
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ex') o (u, U-j ct)-} - o 


- 0 

U x (a/ 4 >) * 0 

- 0 

U^a/b) - i 

u>(e/k,) ~ 0 

* 0 

u> y (s/b)'-o 

«^(a/b) O 

Yields fundamental solutior 

1 ** fc), ^5 <*)j “K (*) 

(V) OCu ( \J' I u>) - O 

u(a/b) - 0 

Ujf^/b) - 0 


!>„ ( a/fe) 0 

u>(a/b) - 0 

<4/3 /fa) -- 1 

(a/b) - 0 

4 ,/s/b) - D 

Yields fundamental solution U t (x) J 

( 4 } 0 'C.U, V} U>^ 0 

u»(a/ b'i 0 

0 

*5 

\ 

ID 

Nw/ 

X 

3 

u(a/b) 

t> x (a/b) 0 

u 3 -(.a/b} " O 

^ xx (a/bi - 0 

u^Ca/fc^ '• o 

u3- (_a/b) ' 1 
xxxr- 


( 6 . 4 . 2 ) 


( 6 . 4 . 3 ) 


( 6 . 4 . 4 ) 


Yields fundamental solution U4 (x)^ TJ^ , (x) 
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In the above described initial value problems the linear 
operator 0(.U,\>,u3-} - o represents the differential equations 
from Eqs. (6.2.31), (6.2.32), and (6.2.33) or Eqs. (6.3.19), 

(6.3.20), and (6.3.21). This operator is also a function of 
the unknown parameter containing the natural frequency of 
vibration. 

If we denote a fundamental solution vector by 



C>7(v) ) 

usrtx)) 


< - 


then by superposition we obtain a general solution 


$ (*) 


T! A; f-Cx) 

i - I 


of the differential equations 


(6.4.5) 


(6.4.6) 


OC&) - O 


(6.4.7) 


This general solution has four arbitrary constants to be 
determined so that the free edge conditions are satisfied. 
The free edge conditions were given by Eqs. (6.2.35) for the 
case of meridional, tangential, and normal deflections and 
by these same equations with a suitable modification of 
the definition of *S r and in the case of radial, 

tangential, and axial deflections. In either case, the 
boundary conditions can be written in operator form 
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o 


( 6 . 4 . 8 ) 

Where B is the linear operator appropriate to the particular 
boundary conditions used. The boundary conditions at the 
free edge are seen to lead to four simultaneous homogeneous 
algebraic equations for the A^ 


.2 A. &(£•)„, - o 


I «• J 


( 6 . 4 . 9 ) 


To emphasize the dependence of the fundamental solutions on 
the frequency parameter (£) we write the algebraic equations 


X a- a 



O 


( 6 . 4 . 10 ) 


These equations will have a non-trivial solution only if the 
determinant of the coefficients of the A^ vanishes. This 
determinant is seen to be a function of the frequency 
parameter (&> so that we will have for a non-trivial solution 
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(6.4.11) 

where is the determinant of the coefficients of the 

set of algebraic equations in Eqs. (6.4.10). 

The eigenvalue problem is seen to be that of finding 
values of for which Eq. (6.4.11) is satisfied. The 

method of approach for solving the eigenvalue equation is 
to evaluate the determinant for a sequence of values of 
(«,) until a sign change in noted. At this point 

a suitable iteration scheme is employed to obtain an 
accurate estimate of the eigenvalue. This scheme is the 
same as the one used in the disk vibration problems but 
is complicated by the fact that each evaluation of the 
determinant requires four numerical integrations of the 
governing differential equations. 

Since we are dealing here with shallow shells we do not 
anticipate the requirement that certain fundamental 
solutions be suppressed to prevent unbounded growth as 
discussed by Zarghamee and Robinson. However, as these 
authors also point out, this feature may have to be available 
if extremely thin shells are to be satisfactorily treated, 
and even this capability may be insufficient if the limit- 
ing case of a true membrane shell is approached. 
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We currently have available the sub-programs required 
for the construction of a digital computer program to carry 
out the above described computations. They include a program 
to compute and tabulate the equilibrium stress and dis- 
placement distribution, a numerical integration scheme for 
determining the fundamental solutions for specified value 
of C%o) , a determinant evaluation procedure for obtaining 
dkCOto) t and an iteration procedure for obtaining the zeroes 
of . In addition, these sub-programs could be used 

to compute the mode shapes for a given frequency if required. 

We have found the available IBM 1620 Computer Facility 
to be too limited to handle problems of this size. While 
most of the sub-programs have been checked out on this 
facility no results have been obtained for the complete 
problem. It is hoped that within the near future these 
computations can be completed at a larger computation 
facility. 
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Figure 1. Geometry of the 


Shallow Spherical Shell. 


L84- 


X 


Figure 2. Coordinate System in which Displacements are 
Measured in the Meridional, Tangential, and Normal 
Directions . 
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Figure 3. Coordinate System in which Displacements are 
Measured Radially, Tangentially, and Axially. 
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Figure 4. Geometry of the Inertia Load. 
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(a) HUB FULLY CLAMPED AT RADIUS r= a 

I 



(b) ANNULAR DISK WITH INNER RADIUS r=Q LOOSELY CLAMPED BY HUB OR RADIUS r=C 


I _ I 



(c) ANNULAR DISK WITH INNER RADIUS r=0. CLAMPED BY HUB OF r=C SO THAT FULL 
CLAMPING EXISTS UP TO r= fir AND PARTIAL CLAMPING. EXISTS FOR £2=rs=C 


Figure 5. Hub Configurations 
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Dimensionless Frequency Parameter 



2 

Figure 8. Variation of the Frequency Parameter, y , as a 
Function of the Hub to Disk Radius Ratio, e, for Various 
Values of the Annulus Radius Ratio, 6, for the Symmetric 
Vibration Case, s=0 , and for the Cases of Zero and One 
Nodal Circles, n=0 and 1* • 
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Figure 12. Variation of the Frequency Parameter, \i , as a 
Function of the Hub to Disk Radius Ratio, e, for Various 
Values of the Annulus Radius Ratio, <5, for the Case of 
One Nodal Diameter, s=l, and for the Case of Two Nodal 
Circles, n=2. 
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Figure 13. Variation of the Frequency Parameter, y , as 
a Function of the Hub to Disk Radius Ratio, e, for Various 
Values of the Annulus Radius Ratio, <5, for the Case of One 
Nodal Diameter, s=l, and for the Case of Three Nodal Circles 
n=3 . 
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Figure 14. Variation of the Frequency Parameter , y , as a 
Function of the Hub to Disk Radius Ratio, e, for Various 
Values of the Annulus Radius Ratio, 6, for the Case of 
Two Nodal Diameters, s=2 , and for the Cases of Zero of One 
Nodal Circles, n=0 and 1. 
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Figure 15. Variation of the Frequency Parameter, y , as a 

Function of the Hub to Disk Radius Ratio, e , for Various 

Values of the Annulus Radius Ratio, 6, for the Case of 

Two Nodal Diameters, s=2 , and for the Case of Two Nodal 

Circles, n=2 . 
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Figure 18.. Variation of the Frequency Parameter, u , as a 
Function of the Hub to Disk Radius Ratio, e, for Various 
Values of the Annulus Radius Ratio, <5, for the Case of . 
Three Nodal Diameters, s=3, and for the Case of Two Nodal 
Circles, n=2. 
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Figure 19. Variation of the Frequency Parameter, v , as a 
Function of the Hub to Disk Radius Ratio> e , for Various 
Values of the Annulus Radius Ratio, 6 , for the Case of Three 
Nodal Diameters, s=3, and for the Case of Three Nodal 
Circles ,n=3 . 
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Figure 21. Tangential Direct Stress Resultant for Freely 
Spinning Shell. Linear Theory with y =0.20. 
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Figure 22. Radial Bending Stress Resultant for Freely 
Spinning Shell. Linear Theory with )) =0.20. 
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Figure 23. Tangential Bending Stress Resultant for Freely 
Spinning Shell. Linear Theory with V =0.20. 
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Figure 25. Radial Direct Stress Resultant for Spinning 
Shell with Fully Clamped Hub. Linear and Nonlinear Theory 
with X =0.058 and V =0.20. 
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Figure 26. Radial Direct Stress Resultant for Spinning 
Shell with Fully Clamped Hub. Linear and Nonlinear Theory 
with X =1.0 and V =0.20. 
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Figure 27. Radial Direct Stress Resultant for Spinning Shell 
with Fully Clamped Hub. Linear and Nonlinear Theory with 
A =3.0 and V =0.20. 
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Figure 28. Radial Direct Stress Resultant for Spinning Shell 
with Fully Clamped Hub. Linear and Nonlinear Theory with 
\ =7.0 and V =0.20. 
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Figure 29. Tangential Direct Stress Resultant for Spinning 
Shell with Fully Clamped Hub. Linear and Nonlinear Theory 
with A =0.058 and V =0.20. 
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Figure 30. Tangential Direct Stress Resultant for Spinning 
Shell with Fully Clamped Hub. Linear and Nonlinear Theory 
with A =1.0 and V* =0.20. 
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Figure 32. Tangential Direct Stress Resultant for Spinning 
Shell with Fully Clamped Hub. Linear and Nonlinear Theory 
with A =7.0 and V =0.20. 
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Figure 34. Radial Bending Stress Resultant for Spinning 
Shell with Fully Clamped Hub. Linear and Nonlinear Theory 
with A =1.0 and V =0.20. 
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Figure 36. Radial Bending Stress Resultant for Spinning 
Shell with Fully Clamped Hub. Linear and Nonlinear Theory 
with A =7.0 and V =0.20. 
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Figure 37. Tangential Bending Stress Resultant for Spinning 
Shell with Fully Clamped Hub. Linear and Nonlinear Theory 
with ^ =0.058 and V =0.20. 
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Figure 38. Tangential Bending Stress Resultant for Spinning 
Shell with Fully Clamped Hub. Linear and Nonlinear Theory 
with A =1.0 and V =0.20. 
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Figure 41. Slope of Deflection Curve for Spinning Shell 
with Fully Clamped Hub. Nonlinear Theory with A =1.0 
and V =0.20. 
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Figure 43. Slope of Deflection Curve for Spinning Shell 
with Fully Clamped Hub. Nonlinear Theory with X =7.0 
and V =0.20. 
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